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Abstract Heuristics indicate that point processes exhibiting clus- 
tering of points have larger critical radius for the percolation 
of their continuum percolation models than spatially homogeneous 
point processes. It has already been shown, and we reaffirm it in this 
paper, that the dcx ordering of point processes is suitable to com- 
pare their clustering tendencies. Hence, it was tempting to conjecture 
that Tc is increasing in dcx order. Some numerical evidences support 
this conjecture for a special class of point processes, called perturbed 
lattices, which are "toy models" for determinantal and permanental 
point processes. However the conjecture is not true in full generality, 
since one can construct a Cox point process with degenerate critical 
radius = 0, that is dcx larger than a given homogeneous Poisson 
point process, for which this radius is known to be strictly positive. 
Nevertheless, the aforementioned monotonicity in dcx order can be 
proved, for a nonstandard critical radius Tc (larger than rc), related 
to the Peierls argument. Moreover, we show the reverse monotonic- 
ity for another nonstandard critical radius (smaller than r^). This 
gives uniform lower and upper bounds on Tc for all point processes 
dcx smaller than some given point process. Moreover, we show that 
point processes dcx smaller than a homogeneous Poisson point pro- 
cess admit uniformly non-degenerate lower and upper bounds on Tc. 
In fact, all the above results hold under weaker assumptions on the 
ordering of void probabilities or factorial moment measures only. Ex- 
amples of point processes comparable to Poisson point processes in 
this weaker sense include determinantal and permanental point pro- 
cesses with trace-class integral kernels. More generally, we show that 
point processes dcx smaller than a homogeneous Poisson point pro- 
cesses exhibit phase transitions in certain percolation models based 
on the level-sets of additive shot-noise fields of these point process. 
Examples of such models are fc-percolation and SINR-percolation. 
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1. Introduction. 

Heuristic. Consider a point process <I> in the d-dimensional Euclidean space 

a given "radius" r > 0, let us join by an edge any two points of 
which are at most at a distance of 2r from each other. Existence of an infinite 
component in the resulting graph is called percolation of the continuum 
model based on Clustering of <I> roughly means that the points of ^ lie 
in clusters (groups) with the clusters being well spaced out. When trying to 
find the minimal r for which the continuum model based on $ percolates, 
we observe that points lying in the same cluster of $ will be connected by 
edges for some smaller r but points in different clusters need a relatively 
higher r for having edges between them. Moreover, percolation cannot be 
achieved without edges between some points of different clusters. It seems to 
be evident that spreading points from clusters of $ "more homogeneously" in 
the space would result in a decrease of the radius r for which the percolation 
takes place. This is a heuristic explanation why clustering in a point process 
$ should increase the critical radius rc = rc{^) for the percolation of the 
continuum percolation model on called also the Gilbert's disk graph or 
the Boolean model with fixed spherical grains. 

Clustering and dcx order. To make a formal conjecture out of the above 
heuristic, one needs to adopt a tool to compare clustering properties of 
different point processes. In this regard, we use directionally convex {dcx) 
order in this article. The dcx order of random vectors is an integral order 
generated by twice differentiable functions with all their second order partial 
derivatives being non- negative.^ Its extension to point processes consists in 
comparison of vectors of number of points in every possible finite collection 
of bounded Borel subsets of the space. Our choice has its roots in Bla- 
szczyszyn and Yogeshwaran [2009], where one shows various results as well 
as examples indicating that the dcx order on point processes implies ordering 
of several well-known clustering characteristics in spatial statistics such as 
Ripley's K-function and second moment densities. Namely, a point process 
that is larger in the dcx order exhibits more clustering, while having the 
same mean number of points in any given set. 

^We remark here that dcx order was initially developed for random vectors (Meester 
and Shanthikumar [1993, 1999], Shaked and Shanthikumar [1990]) partially in conjunction 
with Ross-type conjectures, which predicted that queues with a variable input perform 
worse (cf Miyoshi and Rolski [2004]). Much earlier to these works, a comparative study of 
queues using the supermodular order and motivated by neuron-firing models can be found 
in Buffer [1984]. 
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Conjecture. The above discussion tempts one to conjecture that Tc is in- 
creasing with respect to the dcx ordering of the underlying point processes; 
i.e., $1 <dcx ^2 impUes rc($i) < rc($2)- Since the critical radius rc for 
the percolation does not seem to have any explicit representation as a dcx 
function of vectors of number of points in a finite collection of bounded 
Borel subsets, we are not able to compare straightforwardly Vc of dcx or- 
dered point processes. At the same time, numerical evidences gathered for 
a certain class of point processes, called perturbed lattice point processes, 
were supportive of this conjecture. But as it turns out, the conjecture is not 
true in full generality and we will present a counter-example. 

Counterexample. More specifically, for a given Poisson point process, we 
will construct a dcx larger than it Poisson-Poisson cluster point process (a 
special case of doubly-stochastic Poisson, called also Cox, point process), 
whose critical radius is null (hence smaller than that of the given Poisson 
point process, which is known to be positive). In this Poisson-Poisson cluster 
point process points concentrate on some carefully chosen larger-scale struc- 
ture, which itself has good percolation properties. In this case, the points 
concentrate in clusters, however we cannot say that clusters are well spaced 
out. Hence, this example does not contradict our initial heuristic explanation 
of why clustering in a point process should increase the critical radius for the 
percolation. It reveals rather, that dcx ordering, while being able to com- 
pare the clustering tendency of point processes, is not capable of comparing 
macroscopic structures of clusters. 

Comparison results. What we are able to do with dcx order is to com- 
pare some nonstandard critical radii fc and related, respectively, to the 
finiteness of the expected number of void circuits around the origin and 
asymptotic of the expected number of long occupied paths from the origin 
in suitable discrete approximations of the continuum model. These new crit- 
ical radii sandwich the "true" one, r^{^) < rd^) < fd^). We show that if 
^1 ^dcx ^2, then as suggested by the heuristic, rc($i) < rc($2)- However, 
we obtain the reversed inequality for the other nonstandard critical radius: 
Lci^i) ^ Lc{^2)- This reversed inequality can be also explained heuristically, 
by noting that whenever there is at least one path of some given length in 
a point process that clusters more, there might be actually so many such 
paths, that the inequality for the expected numbers of paths arc reversed. 
However, as we have mentioned above, the fact is that there are examples of 
dcx ordered point processes, for which the inequality is reversed for the crit- 
ical radius Vc itself, namely examples of $i <dcx ^2 with r^^i) > t'c(^*2)- 
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Non-trivial phase transitions. A more positive interpretation of the afore- 
mentioned results is as follows. Combining the two opposite inequalities, we 
obtain as a corollary that if $i <dcx ^2 then the usual critical radius of the 
dcx smaller point process can be sandwiched between the two nonstandard 
critical radii of the larger one, 

rc($2) < r^i^i) < rc($i) < rc($i) < rc($2). 

In more loose terms, the above result can be rephrased as follows: if a point 
process $2 exhibits the following double phase transition < r^{^2) < 
^c(^2) < 00, stronger than the usual one (0 < rc{^2) < 00), then its two 
nonstandard critical radii act as non-degenerate bounds on the critical radius 
rc($), uniformly valid over all $ <dcx ^2, thus in particular implying the 
usual phase transition < rc($) < 00 for all such 

One would naturally like to verify the aforementioned double phase transi- 
tions for the homogeneous Poisson point process. A direct verification eludes 
us at the moment. However, using a slightly modified approach, we are able 
to obtain the same result (uniform bounds on the critical radius Vc) for all 
homogeneous point processes that are dcx smaller than the Poisson point 
process of a given intensity; we call them homogeneous suh-Poisson point 
processes. 

In fact, all the above results regarding comparison of r^, and fc, as well 
as the non-trivial phase-transition for homogeneous sub- Poisson points pro- 
cesses, can be proved under a weaker assumption than the dcx assumption. 
Namely, it is enough to assume that the void probabilities and factorial mo- 
ment measures of the respective point processes are ordered. However, dxc 
ordering needs to be assumed to show non-trivial phase transitions for other 
models based on level-sets of additive or extremal shot-noise fields. 

Examples. In this article, we also provide new examples of point processes 
comparable in dcx order. In particular, families of point processes monotone 
in dcx order are constructed using some random perturbations of lattices. 
Numerical experiments performed for these processes suggest monotonicity 
of Tc with respect to dcx order within this class. 

Recently, determinantal and permanental point processes have been at- 
tracting a lot of attention. They are known as examples of point processes 
that cluster, respectively, less and more than the Poisson point process of 
the same intensity. Assuming trace-class integral kernels, it is relatively easy 
to show that these point processes, have void probabilities and moment 
measures comparable to those of the Poisson point process. Thus, in par- 
ticular, determinantal point processes admit uniform bounds on the critical 
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radius Tc discussed above. Moreover, we show that determinantal and per- 
manental point processes are dcx comparable to the Poisson point process 
on mutually disjoint simultaneously observable sets. Our proof of this lat- 
ter result is based on the known representation of the distribution of these 
point processes on such sets, as well as on Lemma A.l regarding dcx ordering 
of random vectors, which we believe is new, albeit similar to [Meester and 
Shanthikumar, 1993, Lemma 2.17]. Interestingly enough, our perturbed lat- 
tice examples admit a very similar representation, but this time on all sets, 
and it is the same Lemma A.l that allows us to prove their dcx ordering. 

Related work. Let us now make some remarks on other comparison studies 
in continuum percolation. Most of the results regard comparison of differ- 
ent models driven by the same (usually Poisson) point process. In Jonasson 
[2001], it was shown that the critical intensity for percolation of the Boolean 
model on the plane is minimized when the shape of the typical grain is a tri- 
angle and maximized when it is a centrally symmetric set. Similar result was 
proved in Roy and Tanemura [2002] using more probabilistic arguments for 
the case when the shapes are taken over the set of all polygons. This idea 
was also used for comparison of percolation models with different shapes 
in three dimensions. It is known for many discrete graphs that bond per- 
colation is strictly easier than site percolation. A similar result as well as 
strict inequalities for spread-out connections in random connection model 
has been proved in Franceschetti et al. [2005, 2010]. The underlying point 
process in all the above studies has been a Poisson point process and it is 
here that our study differs from them. 

Critical radius of the continuum percolation model on the hexagonal lat- 
tice perturbed by the Brownian motion is studied in a recent pre-print Ben- 
jamin! and Stauffer [2011]. This is an example of our perturbed lattice and 
as such it is a dcx sub-Poisson point process.^ It is shown that for a short 
enough time of the evolution of the Brownian motion the critical radius is 
not larger than that of the non-perturbed lattice, and for large times, the 
critical radius asymptotically converges to that of the homogeneous Poisson 
point process. This latter result is shown by some coupling in the sense of 
set inclusion of point processes. Many other inequalities in percolation the- 
ory depend on such coupling arguments (cf. e.g. Liggett et al. [1997]), which 
for obvious reasons are not suited to comparison of point processes with the 
same mean measures. For studies of this type, convex orders from the theory 

^More precisely, at any time t of the evolution of the Brownian motion, it is dcx smaller 
than a non-homogeneous Poisson point process of some intensity which depends on t, and 
converges to the homogeneous one for t — >■ oo. 
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of stochastic ordering turn out to be quite useful. Our general goal in this 
article is to show the utility of these tools for comparison of properties of 
continuum percolation models. 

Summary of results. Let us be more specific regarding the contributions of 
this paper and quickly summarize our results. 

• Point processes smaller in dcx order have smaller void probabilities as 
well as imply smaller void probabilities of their Boolean models (Propo- 
sitions 3.3 and 3.4). 

• For $1 <dcx ^2, we have rc($i) < fc{^2) and r^{^i) > r^{^2), where 

and fc are certain nonstandard critical radii bounding the (standard) 
critical radius Tc, r^, < Vc <rc, for the percolation of the Boolean model 
with spherical grains (Corollary 4.2 and Proposition 4.5). 

• For any stationary point process having smaller void probabilities and 
factorial moment measures than the Poisson point process of the same 
intensity, < Tc < oo. Moreover, the critical radii rc are uniformly 
bounded away of and of oo over all such point processes of a given 
intensity (Propositions 5.8 and 5.9). 

• The above result is extended to the k— percolation of the Boolean model 
(percolation of the subset covered by at least k of its grains) under the 
assumption that the point process is sub-Poisson (Proposition 5.3). 

• We prove the existence of the percolation regime in certain Signal- 
to-interference- and-noise ratio coverage models driven by homogeneous 
sub-Poisson point processes (Propositions 5.8 and 5.9). These are ex- 
amples of germ-grain models with grains jointly dependent on certain 
shot-noise fields; their percolation has been previously studied assuming 
a Poisson point process of germs. 

• We prove that perturbed lattices with convexly ordered point replication 
kernels are dcx ordered (Proposition 6.1). This provides examples of 
dcx -monotone families of point processes comparable to a given Poisson 
process (Example 6.4). 

• We give an example of a Cox point process with degenerate critical ra- 
dius, rc = 0, which is dcx larger than a given homogeneous Poisson 
process (Example 6.6). 

• We show that determinantal and permanental point processes with trace- 
class integral kernels have void probabilities and factorial moment mea- 
sures ordered with respect to Poisson point processes of the same mean 
measure. Moreover, we show that they are dcx-ordered with respect to 
such Poisson point processes on mutually disjoint simultaneously ob- 
servable sets. (Examples 6.9 and 6.10). 
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Paper organization. The necessary notions and notations are introduced in 
Section 2. In Section 3, we argue that dcx ordering is suitable for the com- 
parison of clustering properties of point processes. In Section 4, we state and 
prove our comparison results for two nonstandard critical radii r^ and fc- 
In Section 5, we use generic exponential estimates to obtain non-degenerate 
bounds for the percolation of various percolation models driven by homo- 
geneous sub-Poisson point processes. New examples of dcx ordered point 
processes, including the counter-example to the conjecture, are provided in 
Section 6. Lemma A.l used for showing dcx ordering of perturbed lattices 
and determinantal and permanental point processes is proved in the Ap- 
pendix. 

2. Notions and notation. 

2.1. Point processes. Let B"^ be the Borel a-algebra and be the a- ring 
of bounded (i.e., of compact closure) Borel subsets (bBs) in the d-dimensional 
Euclidean space R'^. Let M"^ = M(M'^) be the space of non-negative Radon 
(i.e., finite on bounded sets) measures on M'^. The Borel u-algebra Ai"^ is 
generated by the mappings /i i— )• fJ-{B) for all B bBs. A random measure A is a 
random element in (M"^, M'^) i.e, a measurable map from a probability space 
(f], J", P) to (M"^, A^'^). We shah call a random measure <I> a point process 
if <^ G N'^, the subset of counting measures in M'^. The corresponding a- 
algebra is denoted by Af^. Further, we shall say that a point process (pp) $ 
is simple if a.s. $({x}) < 1 for all x G M.'^. Throughout, we shall use A for an 
arbitrary random measure and ^> for a pp. This is the standard framweork 
for random measures (see Kallenberg [1983]). 

As always, a pp or a random measure on M'^ is said to be stationary if its 
distribution is invariant with respect to translation by vectors in W^. 

2.2. Boolean model. Given (the distribution of) a random closed set G 
(racs, see Matheron [1975] for the standard measure-theoretic framework), 
and a pp <I>, by a Boolean model with the pp of germs $ and the typical 
grain G, we call the random set C(<I>, G) = \Jx^(^^{^i + Gi}, where x + A = 
{x + A : a ^ A}, a G M*^, A C M"^ and given <I>, {Gj} is a sequence of i.i.d. 
racs, distributed as G. 

A commonly made technical assumption about the distribution of $ and 
G is that for any compact K (ZW^, the expected number of germs Xi such 
that [Xi + Gi) n K 7^ is finite. This assumption, called for short "local 
finiteness of the Boolean model" guarantees in particular that C($,G) is a 
racs. All the Boolean models considered in this article will be assumed to 
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have the local finiteness property. 

Denote by Bx{r), the ball of radius r centered at x. We shall call G a 
fixed grain if there exists a closed set B such that G = B a.s.. In the case 
B = Bo{r), where O is the origin of Mf^ and r > constant, we shall denote 
the Boolean model by C{^, r). Note that C(<I>, r) is a racs and so is C{^, G) 
for any G a.s. bounded. 

2.3. Directionally convex ordering. Let us quickly introduce the theory 
of directionally convex ordering. We refer the reader to [Miiller and Stoyan, 
2002, Section 3.12] for a more detailed introduction. 

The order < on M*^ will denote the component- wise partial order, i.e., 
(xi, . . . , Xk) < (yi, . . . , Vk) if Xi < yi for every i. For p, q,x,y G M'^ we shall 
abbreviate p < min{x, y} to p < {x, y} and max{x, y} < q to {x, y} < q. 

One introduces the following families of Lebesgue-measurable functions on 
M'^: a function / : M'^ — t- M is said to be directionally convex (dcx) if for every 
X, y,P,q G IS'^ such that p < {x, y} < q and x + y = p + q, we have that 

fix) + fiy)<f{p) + fiq). 

A function / is said to be directionally concave (dcv) if the inequality in 
the last equation is reversed. A function is said to be supermodular (sm) 
if the above inequality holds for all x, y G M'^' with p = min{a;, y} and 
q = max{x,?/}. A Lebesgue-measurable function is dcx if and only if it is 
sm and coordinate- wise convex. A convex function on M will be denoted by 
ex. We abbreviate increasing and dcx by idcx., decreasing and dcx by ddcx 
and similarly for dcv functions. 

Let 5^ denote some class of Lebesgue-measurable functions from M*^ to M 
with the dimension k being understood from the context. In the remaining 
part of the article, we will mainly consider ^ to be one among the class of 
dcx , idcx , idcv , dcv , ddcv , ddcx functions. Unless mentioned, when we state 
E{f{X)) for / G 5^ and X a random vector, we assume that the expectation 
exists. Suppose X and Y are real-valued random vectors of the same dimen- 
sion. Then X is said to be less than Y in ^ order if E{f{X)) < E{f{Y)) for 
all / G 5^ such that both the expectations are finite. We shall denote it as 
X <5 Y. This property clearly regards only the distributions of X and Y, 
and hence sometimes we will say that the law of X is less in ^ order than 
that of Y. 

Suppose {^(s)}^;^^ and {y(s)}sG5 are real-valued random fields, where 
S is an arbitrary index set. We say that {X{s)} <^ {Y{.s)} if for every k > 1 
and si, . . . ,Sk £ S, {X{si), X{sk)) <s iY{si), Y{sk)). 
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A random measure A on M"' can be viewed as the random field {A{B)}^^^d. 
With the aforementioned definition of dcx ordering for random fields, for two 
random measures on M'^, one says that Ai(-) <dcx ^2(")) ^-^y -^Ij • • • > 
bBs in M'^, 

(1) (Ai(Si),...,Ai(Sfc)) <dc, (A2(5i),...,A2(5fc)). 

The definition is similar for other orders, i.e., when 5 is the class of idcx, idcv, 
ddcx or ddcv functions. It was shown in Blaszczyszyn and Yogeshwaran 
[2009] that it is enough to verify the above condition for Bi mutually disjoint. 

In order to avoid technical difficulties, we will consider here only ran- 
dom measures (and pp) whose mean measures E(A(-)) are Radon (finite on 
bounded sets). For such random measures, dcx order is a transitive order ^. 
Note also that Ai(-) <dcx ^2(") implies the equality of their mean measures: 
E(Ai(-)) = E(A2(-)). For more details on dcx ordering of pp and random 
measures, see Blaszczyszyn and Yogeshwaran [2009]. 

3. Clustering and dcx ordering of point processes. In this section, 
we will present some basic results on dcx ordering of pp that will allow us 
to see this order as a tool to compare clustering of pp. It will be shown 
that pp smaller in dcx order clusters less. These results involve comparison 
of some usual statistical descriptors of spatial (non-)homogeneity as well as 
void probabilities. Examples of comparable pp will be provided in Section 3.4 
and more extensively in Section 6. 

3.1. Statistical descriptors of spatial (non-) homogeneity. Looking at Fig- 
ure 1, it is intuitively obvious that some pp cluster less than others. However, 
to the best of our knowledge, there has been no mathematical formalization 
of such a statement based on the theory of stochastic ordering. In what fol- 
lows, we present a few reasons to believe that pp that are smaller in dcx 
order exhibit less clustering. 

Roughly speaking, a set of points in is "spatially homogeneous" if ap- 
proximately the same number of points occurs in any circular region of a 
given area. A set of points "clusters" if it lacks spatial homogeneity. There 
exist some statistical descriptors of spatial homogeneity or clustering of pp. 
We will show that dcx ordering is consistent with these statistical descrip- 
tors. We restrict ourselves to the stationary setting. 

^Due to the fact that each dcx function can be monotonicaUy approximated by dcx 
functions fi{-) which satisfy fi{x) = O(||a;||oo) at infinity, where ||a;||oo is the Loo norm on 
the Euclidean space; cf. [Miiller and Stoyan, 2002, Theorem 3.12.7]. 
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simple perturbed lattice 



Poisson point process 



Cox point process 



Figure 1. From left to right : patterns sampled from a simple i.i.d. perturbed lattice, 
Poisson pp and a doubly stochastic Poisson (Cox) pp, all with the same mean number of 
points per unit of surface. For more details see Example 6.4 in Section 6.1. 

One of the most popular functions for the statistical analysis of spatial 
homogeneity is the Ripley's K function defined for stationary pp (cf Stoyan 
et al. [1995]). Assume that $ is a stationary pp on M"^ with finite intensity 
A = E($([0,1]°')). Then 

K{r):=^p[ mBxAr))-l)). 

where ||-B|| denotes the Lebesgue measure of a bBs B. Due to stationarity, 
the definition does not depend on the choice of B. The following observation 
was made in Blaszczyszyn and Yogeshwaran [2009]. 

Proposition 3.1. Consider two stationary pp ^j, j = 1,2, with the 
same finite intensity and denote by Kj{r) their Ripley's K functions. If 
<dcx ^2 thenKii-) < K2{-). 

Another useful characteristic for measuring clustering effect in point pro- 
cesses is the pair correlation function. It is related to the probability of 
finding the center of a particle at a given distance from the center of an- 
other particle and can be defined on M?'^ as g{x,y) := p{?) (^JfpW {y) ' '^^^^'^^ 
p^^^ is the kth joint intensity; i.e., a function (if it exists) p^^'^ : (R'^)^ — )• 
[0,00) such that for any k mutually disjoint subsets Bi, . . . , Bk, we have 
that Enf=i ^{Bi) = fB^x-xBk P^'^^i^^' ■ ■ ■,Xk)dxi . ..dxk (see Stoyan et al. 
[1995]). 

We now make a digression to joint intensities, which will be used later on. 
Recall that the measure a''{-) defined by a^(i?i x . . . x 5^) = EH^Li ^{Bi) 
for all (not necessarily disjoint) bBs Bi [i = 1,...,A;) is called kth order 
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moment measure of For simple pp, the truncation of the measure a^{-) 
to the subset {{xi, . . . , Xk) S (M'^)^ : Xi 7^ xj, for i ^ j} is equal to the k th 
order factorial moment measure a^^\-). Hence p^^'^ (if it exists) is the density 
of a^'^^(-) for simple pp. Consequently, the joint intensities characterize the 
distribution of a simple pp. The above facts remain true even when the 
densities are considered with respect to HiLi IJ-idxk) for an arbitrary 
Radon measure /i on M"^. We will assume this generality in Examples 6.9 
and 6.10. 

The following result was also proved in Blaszczyszyn and Yogeshwaran 
[2009]; 

Proposition 3.2. Let $1, $2 he two pp on with a-finite k th moment 
measures a^, j = 1,2. If <idcx $2 then a^(^) < a^{A) for all bBs 

A C (M'^)*^. Consequently pf^\xi, . . . ,Xk) < p[''\xi, . . . ,Xk) (whenever they 
exist) for Lehesgue a.e. {xi . . . ,Xk)- Moreover, if^i <dcx ^2, then gi{x,y) < 
92{x,y) for Lehesgue a.e. {x,y). 

Recall that for a (stationary) Poisson pp of intensity A with < A < 
00, we have that K{r) = 1 and p{x, y) = 1. Poisson pp is often considered as 
a reference model of spatial homogeneity to which other pp are compared. 
In Sections 3.3 and 6, we will consider pp smaller and larger than ^\ in dcx 
ordering. 

3.2. Void probabilities . One of the fundamental observations to our ap- 
proach is that dcx ordering allows us to compare the void probabilities of 
pp y[B) = P ($ n S = 0) for all bBs B. 

Proposition 3.3. Denote by i^i{-),V2{-) the void probabilities of pp <J>i 
and $2 on respectively. If <ddcx ^2 then vi{B) < V2{B) for all bBs 
B C W^. 

Proof. This follows directly from the definition of dcx ordering of pp, ex- 
pressing i^j{B) = E(/($j(i?))), j = 1, 2, with the function f{x) = max(0, 1 — 
x) that is decreasing and convex (so ddcx in one dimension). □ 

The latter result can be interpreted as follows: for two pp that are ddcx 
ordered, the smaller one has less chance to create a particular hole (absence 
of points in a given region). In the case of dcx ordered pp this holds true 



*(j-finiteness condition is missing there; see [Yogeshwaran, 2010, Prop. 4.2.4] for the 
correction. 

11 



even if both pp have the same expected number of points in any given set. 
The above result can be strengthened to the comparison of void probabih- 
ties of Boolean models having their germ pp dcx ordered. This observation 
will be the starting point for our investigation of the connections between 
percolation and directionally convex ordering of pp in Section 4. 

Proposition 3.4. Let C{^j,G), j = 1,2 be two Boolean models with pp 
of germs ^j, j = 1,2, respectively, and common distribution of the typical 
grain G. Assume that $j are simple and have locally finite moment measures. 
If^i <dcx $2 then P G) n S = 0) < P (C($2, G) n 5 = 0) for all bBs 

B <ZW^ . Moreover, if G is a fixed compact grain then the same result holds, 
provided i'i{B) < V2{B) for all bBs B C , where Ui is the void probability 
of^i. 

Proof. For xgR'^, F racs and B G B^, define a function h{{x, F),B) = 
l[Bri{x + F) 7^ 0] and an extremal shot-noise field on defined by S^{B) := 
s^P(^X',G')&S>i Ki^p G^j), B), where = {(Xj, G*)} are the pp of germs in- 
dependently marked by the grains defining the corresponding Boolean mod- 
els. By the local finiteness of the two Boolean models, P (C(<I>i, G) n = 0) = 
P {Si{B) < 0) for all bBs B C M'^. By [Blaszczyszyn and Yogeshwaran, 
2009, Proposition 3.2, point 2], $i <dcx ^2 implies l>i <dcx ^2 (^)- More- 
over, by Proposition 4.1 in the aforementioned article, P {Si{B) < 0) < 
P {S2{B) < 0). This completes the proof of the first statement. 

The proof of the result for the Boolean model with a fixed compact grain 
follows immediately from the equivalence: (a; + G) H -B = if and only if 
X ^ B + G, where G = {-y :yeG}. □ 

3.3. Sub- and super-Poisson point processes. We now concentrate on 
comparison of pp to the Poisson pp of same mean measure. To this end, 
we will call a pp dcx sub-Poisson (respectively dcx super-Poisson) if it is 
smaller (larger) in dcx order than the Poisson pp (necessarily of the same 
mean measure). For simplicity, we will just refer to them as sub-Poisson or 
super-Poisson pp omitting the word dcx. Examples of such pp are given in 
Section 6. In particular, a rich class of such pp called the perturbed lat- 
tice pp will be provided later. Also it was observed in Blaszczyszyn and 
Yogeshwaran [2009] that Poisson-Poisson cluster pp, Levy based Cox pp, 
Ising- Poisson cluster pp etc... are super-Poisson pp. 

^The assumption on local-finiteness of mean measures of "l>i is missing in the statement 
of this result in Blaszczyszyn and Yogeshwaran [2009]. Also the idcx / ddcx / idcv ordering 
is mentioned there but not proved. 
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We will also consider some weaker notions of sub- or super-Poisson pp, 
for which only moment measures or void probabilities can be compared. 
Recall from Propositions 3.2 and 3.3 that dcx order implies ordering of 
factorial moment measures and void probabilities z^(-). Recall also that a 
Poisson pp can be characterized as having void probabilities of the form 
= exp{—a{B)). Bearing this in mind, we say that a pp $ is weakly sub- 
Poisson in the sense of void probabilities (z^-weakly sub-Poisson for short) 
if 

(2) P($(B) = 0) <e-^*(^)) 

for all Borel sets B C M'^. When the inequality in (2) is reversed, we will say 
that <I> is weakly super-Poisson in the sense of void probabilities (z^-weakly 
super-Poisson). 

In our study of the percolation properties of pp, it will be enough to be 
able to compare moment measures (off diagonals) of pp. In this regard, we 
say that a pp <^ is weakly sub-Poisson in the sense of moment measures 
(a-weakly sub-Poisson for short) if 

(k \ ^ 

HHbM <llE{^{B,)) 
i=l / i=l 

for all mutually disjoint bBs Bi C M'^. When the inequality in (3) is reversed, 
we will say that ^ is weakly super-Poisson in the sense of moment measures 
{a-weakly super-Poisson). 

Finally, we will say that <1> is weakly sub-Poisson if <I> is a-weakly sub- 
Poisson and z^-weakly sub-Poisson. Similarly, we define weakly super-Poisson 
pp. 

3.4. Simple Examples. We give here some examples of weakly super- 
Poisson pp that follow immediately from earlier results. More detailed ex- 
amples of sub- and super-Poisson pp shall be presented in Section 6. 

It can be seen that any associated pp ^ is a-weakly super-Poisson. In 
fact, it is easy to see using the Jensen's inequality that the associated pp 

®i.e, pp <& satisfying 

(4) Cov . . . , ^{Bk))giHBi), HBk))) > 

for any finite collection of bBs _Bi , . . . , Bk C R"^ and /, g continuous and increasing func- 
tions taking values in [0,1]; Burton and Waymire [1985]. This property is also called 
positive correlations, or the FKG property. 
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have moment measures everywhere (not only off the diagonals) larger than 
those of the corresponding Poisson pp ^. From [Burton and Waymire, 1985, 
Th. 5.2], we know that any Poisson center cluster pp is associated. This is a 
generalization of our perturbation (19) of a Poisson pp $ (cf. Example 6.3) 
having form = X]j^g^{^ + <&x} with <I>x being arbitrary i.i.d. (clus- 

ter) point measures. Other examples of associated pp given in Burton and 
Waymire [1985] are Cox pp with intensity measures being associated. 

It is easy to see by the Jensen's inequality that all Cox pp are v-weakly 
super-Poisson and hence Cox pp with associated intensity measures are 
weakly super-Poisson. 

In Section 6.3, we will show that determinantal and permanental pp are 
weakly sub-Poisson and weakly super-Poisson respectively. Moreover, their 
dcx comparison to Poisson pp is possible on mutually disjoint, simultaneously 
observable sets. 

4. Percolation of Boolean models and the dcx ordering of the 
underlying point processes. By percolation of a Boolean model C = 
C(<I>,G) (see Section 2.2), we refer to the existence of an topologically un- 
bounded connected subset (giant component) of C. We have given heuristic 
arguments in the Introduction as to why more clustering in the pp $ should 
make percolation of the corresponding Boolean model C($,G) less likely. 
Since we have shown in the previous section that dcx order is suitable for 
comparing the clustering tendencies of pp, it is tempting to conjecture that 
pp smaller in dcx order percolate better. Some numerical evidences support- 
ing this conjecture for a certain class of pp, called perturbed lattice pp, will 
be presented in Section 6.1.1. However, as we have also mentioned in the 
Introduction, this conjecture is not true in full generality; a counterexam- 
ple will be presented in Section 6.2. Nevertheless, some general comparison 
results regarding nonstandard critical radii for percolation and dcx order of 
pp shall be proved in this section. 

Define a critical radius for C(<I>,G) as 

(5) rc = rc($, G) := inf |r > : P (c{<i>, G^''^) percolates) > o} ; 

i.e., the smallest r for which C(<I>, G*-^-*) percolates with a positive probability, 
where G^^^ := [JxeG -^^i''^)^ ^"^'^ parallel set of G. 

^Indeed, due to the symmetry of the moment measures it is enough to verify the 

inequahty on rectangles A = B^^ x • ■ • x B^" for any mutually disjoint bBs Bi C K'', any 
ki + . . . + kn — k, ki £ Z, k, ki > 1 with X]"^! ~ ^- Recall that for the Poisson pp of 

mean measure a(-), we have ce''{A) = ^X^jLi (j') (^{Bki)^ ^- 
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We say that the origin percolates in C($,G) if O belongs to an infinite 
component of C($, G). If $ is stationary then it is easy to show that for any 
r > 0, C($,G(''^) percolates with a positive probability if and only if the 
origin percolates with a positive probability. Further, ergodic pp percolate 
with probability that is either zero or one. 

The critical radius rc is similar to the critical radius defined for percolation 
in various continuum percolation models. For example, taking G = {O} we 
obtain the critical radius r for the percolation of the classical continuum 
percolation model (the Boolean model with spherical grains G = Bo{r)). 

In what follows, we will define two other critical radii for the percolation 
of C($, G). They will be an upper and a lower bound of rc respectively. We 
will show that the ordering of void probabilities and moment measures of 
the pp germs $ allows comparison of these two new critical radii. In this 
regard, we introduce some more notations. 

4.1. Auxiliary discrete models. Though wc focus on the percolation of 
Boolean models (continuum percolation models), but as is the wont in the 
subject we shall extensively use discrete percolation models as approxima- 
tions. For r > 0, a; G M*^, define the following subsets of W^. Let := 
(-i>^]^ Q"i^) ■= ^ + Q"^ Qr ■■= {-r,rY and Q,(x) := x + Q,.. We 
will consider the following three families of discrete graphs parametrized by 
n G N. 

• L*'^ = (Z^,E*'^) is the usual close-packed lattice graph scaled down 
by the factor 1/n. It has = ^Z'', where Z is the set of integers, 
as the set of vertices and the set of edges E**^ := {{zi,Zj) G (Z^)^ : 
Q^{z^i)r\Q^{z:j)^^. 

• L*°'(r) = (rZ'^, E*'^(r)) is a similar close-packed graph on the scaled- 
up lattice rZ*^; the edge-set is W^{r) := {{zi,Zj) G {rlJ^f : Qrizi) fl 
QriZj) + 0}. 

• Finally, for a given compact set K C R'', the graph II^{K) = {Zf„ E;'^(K)) 
has vertex-set Z^ and edge-set E*J{K) := {{zi,Zj) : {Q''{zi) + K} n 
{Q''{zj)+K} 7^ 0}, where we recall that A+B = {a+b : a e A,b e B}. 

A path in h*/{K) from the origin O G M*^ to dQm is 7 = (^^i, . . . , Zk+i) C 
Zf^ such that foi I < i < k, Zi e Qm, {zi, z.+i) G e{{K) and O G + 
K}, {Q"(zfe+i) + K}n dQ„, / 0. Let Tl^{K) be the set of all such paths 
from the origin to dQm in I-'*'^(i^). 

A contour in L*'^ is a minimal collection of vertices such that any infinite 
path in L*'^ from the origin has to contain one of these vertices (the mini- 
mality condition implies that the removal of any vertex from the collection 
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will lead to existence of an infinite path from the origin without any inter- 
section with the remaining vertices in the collection). Let r„ be the set of 
all contours around the origin in L*'^. For any subset of points 7 C M"^, in 
particular for paths 7 G T^(i^), Tn, we define = [Jze'y Q^'i^)- 

Throughout the article, we will define several auxiliary site percolation 
models on the above graphs by randomly declaring some of their vertices 
(called also sites) open. As usual, we will say that a given discrete site perco- 
lation model percolates if the corresponding sub-graph consisting of all open 
sites contains an infinite component. In particular, note that the existence 
of only a finite number of paths in r„ (i.e., contours around the origin in 
L*'^) implies the percolation of the corresponding site percolation model on 
L*'^. If this number of paths in r„ has a finite mean (which is thus a suffi- 
cient condition for the percolation) then we will say that the site percolation 
model on L*°' percolates through the Peierls argument (cf. [Grimmett, 1999, 
pp. 17-18]). 

4.2. Ordering of void probabilities and percolation. Define a new critical 
radius 
(6) 

= rc($, G) := infjr > : for aU n > 1, ^ P (c{<i>, G^''^) n = 0) < 00} . 

76r„ 

It might be seen as the critical radius corresponding to the phase transition 
when the discrete model L**^ = (Z^,E*'^), approximating G($,G'^^)) with 
an arbitrary precision, starts percolating through the Peierls argument. As 
a consequence, rc(H) is an upper bound for the actual critical radii, as we 
show now. 

Lemma 4.1. Let G(<I>,G), be a Boolean model with the pp of germs $ 
and G as the distribution of the typical grain. Then, fc(<I>,G) > rc{^,G). 

Proof. The proof as indicated above shall use Peierls argument on a 
suitable discrete approximation of the model on L*'^. Define a random field 
Xr := {Xr{z)]^^j^d such that X,(z) = 1[G(«>, G^''^) n Q"(z) / 0]. By Peierls 
argument for r > r^ for all n, Xr percolates in a.s.. Thus for all n, 
G($, G^^"*"^-*) percolates a.s. and so r + ^ > r^. Hence, r > Tc. □ 

Here is the main observation of this section. It follows directly from the 
definition of fc and Proposition 3.4. 
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Corollary 4.2. Let Cj = C{^j,G), j = 1,2 be two Boolean models 
with pp of germs j = 1,2, respectively, and G as the common distribu- 
tion of the typical grain. Assume that <I>j are simple and have locally finite 
moment measures. <dcx ^2 thenrc{^i,G) <fc{^2,G). Moreover, if G 
is a fixed compact grain then the same result holds provided i'i{B) < 1^2 (-B) 
for all bBs B C , where Vi is the void probability of ^i. 

Remark 4.3. The critical radii Tc = rcC^) and fc = ?^c(S) can be de- 
fined for an arbitrary racs H by replacing by H(^) in (6) and (5), 
respectively. The result of Lemma 4.1 remains true (cf. the proof). More- 
over, if two racs Hi,H2 are ordered with respect to their void probabili- 
ties; i.e., P(Hini^ = 0) < P(H2n/C = 0) for any compact set K, then 



4.3. Comparison of moment measures and percolation. We will define 
now another critical radius for C{^,G) by counting some paths on it. In 
this regard, define an open path 7 from the origin to dQm in G) (recall 
the notation from Section 4.1) as an (ordered) subset {Xi, . . . ,Xk+i) of 
distinct points of <I> such that for all 1 < i < A;, {X.i}i>i C Qm-, {^i+i + 



Gi+i} n {Xi + / and {Xk+i + Gk+i} n OQ^ / , O G {Xi + d}. 



Let Nm{^, G) be the number of such open paths from the origin to dQ 
Define the following new critical radius : 



It corresponds to phase transitions in the Boolean model when the expected 
number of open paths from the origin to the boundary of an arbitrary large 
box becomes positive. It is easy to see that is a lower bound for re- 

Lemma 4.4. Assume a stationary pp <I> of germs and let the typical grain 
G be an almost surely connected set. Then r^{^,G) < rc{^,G). 

Proof. By the stationarity of <I>, for any r > 0, C{^,G^'^^) percolates 
with positive probability if and only if the origin percolates in C{^,G^^'^) 
with positive probability. Moreover, by the connectivity of grains, this latter 
statement is equivalent to the statement that P [Nm{^,G^'^^) > l) > for 
all m > 1. Consequently, since the above sequence of events is decreasing in 
m. 



rcC^i) < rc(H2). 





Now, the proof follows by the Markov's inequality. 



□ 
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Here is the main result of this section. 

Proposition 4.5. Let Cj = C{^j,G), j = 1,2 be two Boolean models 
with the same fixed (deterministic) grain G and simple pp of germs ^j, 
j = 1, 2, with a -finite k th moment measures a^, respectively, for all k > 1. 
// a^{Bi X • • • X Bk) < a^iBi x • • • x B^) for all k > 1 and all mutually 
disjoint Bi, i = 1, . . . , k, then r^{^i, G) > 'Lc{^2-, G). In particular, the above 
inequality is true if ^dcx ^2- 

Note the reversal of the inequality in the above Proposition for the critical 
radii with respect to the result of Corollary 4.2 regarding fc- A similar 
result can also be proved for the critical radius defined with respect to the 
expected number N*, of paths of length n in the Boolean model, as well as 
for the expected number of crossings of a rectangle ^. 

Proof of Proposition 4.5. Define r^ := infjr : E[Nm{^2,G'^''"^)) = 

oo| for any m > 1. (Recall that Nm{^i,G^'^'^) is the number of open paths 
from the origin to dQm in C7($i,GW)). Wewill first prove that the condition 

(8) E(^Nm{1>i, GM)) < E(^Nm{^2, G^^^)) , for aU m > 1 and r / r*^ 

implies r^{^i,G) > r^($2, G). 

Indeed, assume (8) and suppose that Zci^i^ < r^{^2, G). Then choose 
r {?^m}meN such that r^i^i, G) < r < rci^2, G) and note that (8) implies 
< liminf^E(A^^($i,G'W)) < liminf^ E(iV^($2, G^)). This contradicts 
r < r^(^>2, G) and hence r^{<^i,G) > r^{^2, G). 

Note that (8) is true for r > r^. Thus, in order to prove (8), we will work 
with a fixed m > 1 and r < r^ for the remaining part of the proof. 

^Similarly to open paths from the origin to dQm, one can define an open path on the 
germs of $ crossing the rectangle [0, m] x [0, 3m] x . . . x [0, 3m] across the shortest side 
and define yet another critical radius rs as the smallest r for which such a path exists 
with positive probability for an arbitrarily large m. If all grains are path- connected then 
the existence of such a path on the germs of $ is equivalent to the existence of a crossing 
defined as a continuous curve in C("1?,G). This consistency allows us to relate this new 
critical radius rs with the critical intensity defined in [Meester and Roy, 1996, (3.20)]. 
Hence, by [Meester and Roy, 1996, Proposition 3.5], — r^ for Boolean models with 
an homogeneous Poisson pp of germs and spherical grains of bounded radius (possibly 
random). Moreover, replacing the probability of the existence of such crossing by the 
expected number of crossings, we can define a lower bound < in the same way as we 
have defined < rc. 
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Let and ^>2 be the given independent pp of germs and G be the fixed 
(deterministic) grain. Let R be such G C Bo{R), where Bo{R) is the bah 
centered at the origin of radius R. The first observation is that in order to 
study the paths from the origin to dQm on C(<I>i, G^^^), i = 1,2, and on the 
discrete graph L*'^(G*^'^)), n > 1, with < s < r + Vd/2, it is enough to 
consider these models only in Qm' with m' = m + R + r + \fd. 

Let no be the smallest n > 1 such that for all z G Qm' H either 
($1 + $2)(Q"(^)) < 1 or ($1 + ^2)iCt[z)) = 2 and $i and $2 share a 
common atom in Q^{z). The random variable no is finite since both $1 and 
^2 are simple pp. 

Let us now remark the following two useful relations between Boolean 
models C($j, G^'^^), i = 1,2, and their discrete graph approximations h^^{G^^^), 
n > 1. 

• For any n > 1, 

(9) Nmi^i, GW)l(n > no) < iV;^($i, GM)l(n > no) , 

where for a compact set K, N^{^i,K) = E7eTs,(A') H^st 
and (recall from Section 4.1) T^iK) denotes the set of open paths 
from the origin to dQm in the graph L*'^(i^). This follows from the fact 
that on the event {n > no}, for any z G Z^, if ^i{Q^{z)) > then $j 
has at most one point in Q'^{z). Moreover, for any two distinct points 
Xi,X2 € ^i, Xi / X2 , necessarily belonging to two different boxes 
Q^{zi), Q"'{z2), respectively, with zi / Z2, which satisfy {Xi + G(''")}n 
{X2 + / 0, we have that {Q"(zi) + G^^')} n {Q"(z2) + G^'')} / 0. 

• For any n > 1 and s(n) = r + ^/d/{2n), 

(10) iV-($„GW)<iV™(cI>„GW"))). 

This is so because for any X G <I>j n Q^{z) for some 2; G Z^, we have 
{Q"(z) + GW} C {X + G(^("))} and, moreover, in iV;^($i, G^) we are 
counting only the paths of G(<I>j, G*^**^"^^) from O to dQm, which do 
not have edges between points in the same box Q"'{z) for some 2; G Z^. 

The following observation regarding the discrete models is crucial to our 
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proof. 

E(iv;^($i,GM)) = Yl E(iV,($i)) 

7GTS,(GM) 

^ E E/'n^2(Q"(^.))] 

(11) = E(iv;;,($2,GM)) 

where the inequality is due to our assumption on the moment measures. 
Using (9) and Fatou's lemma we get 



Hminf E( Af;^($i,G(''))) > E( hminf iV^($i, G('^))l(n > no)} 
:i2) = EfiV™(<I>i,GM) 



For r < rj^, choose ni such that s(ni) < rj^. Note that for n > ni, again 
by (10) we have 

iV;^(^2,GM) < iV„(cI>2,G(^("))) < iV„(<I>2,G(^("^))) 

and, by the definition of r"^, E[Nm{^2, G^^^"'^^^)) < oo. Hence, using reverse 
Fatou's lemma, (10), (11) and (12), we get that 



E iV„(<Di,GM) < hmsupE iV;^(cI>i,GM 

< limsupEfc($2,G(") 

< E(limsupiV;^($2,G(")) 

< E(limsupiV^($2,G("^"^)) 



E(iV™(cI>2,GM)) , 



where the last equality follows from the fact that (~]^^iC{^2,G^^^"'^^) = 
C{^2,G'^), since Boolean model is a closed set. Thus we have shown (8), 
which concludes the proof. □ 

Remark 4.6. Combining the results of Corollary 4.2 and Proposition 4.5 
regarding percolation of the two Boolean models C($i, G), with a fixed con- 
nected grain G and simple pp of germs of locally finite moment measures, 
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we obtain the following sandwich inequality for the critical radius rd^i, G) 
of the Boolean model C(<I>i, G) 

r,{^2, G) < r,($i, G) < re($i, G) < fc($i, G) < rc(^2, G) , 

provided is stationary and has smaller void probabilities and facto- 
rial moment measures than ^2- Assume moreover that r^{^2^G) > and 
rc{^2,G) < oo. Then the above sandwich inequality implies that C{^i,G) 
exhibits a non-trivial phase transition relative to the percolation radius; i.e., 
< rc(^>i,G) < oo. 

5. Non-trivial phase transition for percolation models on sub- 
Poisson point processes. Following Remark 4.6, with $2 being a homo- 
geneous Poisson pp ^x, one could prove the existence of a non-trivial phase 
transition in the Boolean model driven by a sub-Poisson pp (or weakly sub- 
Poisson pp in the case of fixed grains) <I'i, by showing that r^{^x,G) > 
and rc{^x,G) < 00. Even for a Poisson pp, a direct verification of these 
latter conditions evades us. However, the existence of a nontrivial phase 
transition for the percolation of Boolean model and other models driven by 
sub-Poisson pp can be shown and that shall be the goal of this section. 

We will be particularly interested in percolation models on level-sets of 
additive shot-noise fields. The rough idea is as follows: level- crossing proba- 
bilities for these models can be bounded using Laplace transform of the un- 
derlying pp. For sub-Poisson pp (pp that are dcx smaller than Poisson pp), 
this can further be bounded by the Laplace transform of the corresponding 
Poisson pp, which has a closed-form expression. For 'nice' response func- 
tions of the shot-noise, these expressions are amenable enough to deduce 
the asymptotic bounds on the expected number of closed contours around 
the origin or the expected number of open paths of a given length from 
the origin and thus, using standard arguments, deduce percolation or non- 
percolation of a suitable discrete approximation of the model. In what follows 
we shall carry out this program for /c-percolation in the Boolean model and 
percolation in the SINK model. For the similar study of word percolation 
see [Yogeshwaran, 2010, Section 6.3.3]. 

5.1. Bounds in discrete models. We shall start with a generic bound on 
a discrete model which shall be used to prove bounds in the continuum 
models. Denote by V$(x) := Ylxi^^ ^(^^ (additive) shot-noise field 

generated by a pp <I> and a non- negative response function defined on 
X M'^. Define the corresponding lower and upper level sets of this shot- 
noise field on the lattice rZ"' by Z^(V$, < h) := {z G rZ*^ : V$(z) < h} and 
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^r(^*> > h) := {z € rU^ : V^{z) > h}. We will be interested in percolation 
of Z^(V$, < h) and > /i) understood in the sense of site-percolation 

of the close-packed lattice L*'^(r) (cf Section 4.1). 

Remark 5.1. Recall that the number of contours surrounding the origin 
in L*'*(r) ^ is at most n(3'^ — 2)"~^. Hence, in order to prove percolation 
of a given model using Peierls argument, it is enough to show that the 
corresponding probability of having n distinct sites simultaneously closed is 
smaller than for some < p < {3'^ - 2)-^ for n large enough. Similarly, 
since the number of paths of length n starting from the origin is at most (S*^ — 
1)", in order to disprove percolation of a given model it is enough to show 
that the corresponding probability of having n distinct sites simultaneously 
open is smaller than for some < p < (S'^ — 1)~^ for n large enough. 

The following result allows us to derive the afore-mentioned bounds. We 
restrict ourselves to the stationary case. 

Lemma 5.2. Let ^ be a stationary pp and V$(-), Z^(V$, < h), Z^(V$, > 
h) be as defined above. Let ^\ be the homogeneous Poisson pp with intensity 
A on W^. If ^ <idcx then for any s > 0, 

(13) P(VcE.(zi)>/i,l<i<n)<e-^"''exp|A j [e'^^-^^'^''''^^ - l)dx^ . 
If ^ <ddcx then for any s > 0, 

(14) P{V^{zi)<h,l<i<n)<e'''^exp!.X [ (e""^-! ^(^'^'^ - l)dxl . 

I jR'i J 

Proof. In order to prove the first statement, observe by Chernoff's in- 

contour surrounding the origin in L*'*(r-) is a minimal collections of vertices of 
L*''(r) such that any infinite path on this graph from the origin has to contain one of 
these vertices. 

i°The bounds n(3'' - 2)"-^ and {3'^ - 1)" are not tight; we use them for simplicity of 
exposition. For more about the former bound, refer Balister and BoUobas [2007] , Lebowtiz 
and Mazel [1998]. 
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equality that for any s > 0, 



P{V^{zi)>h,l<i<n) < 

< 
< 



where the third inequality follows from the ordering $ <idcx ^\ ^iid the 
equality by the known representation of the Laplace transform of a functional 
of Poisson pp (cf [Daley and Vere-Jones, 2007, eqn. 9.4.17 p. 60]). 

The proof of the second statement follows along the same lines by not- 
ing that for any random variable X and any a G R, s > 0, P {X < a) = 
P {e-'^ > e"™) < e'"'E{e-''^) . □ 

5.2. k-percolation in Boolean model. By fc-percolation in a Boolean model, 
we understand percolation of the subset of the space covered by at least k 
grains of the Boolean model. The aim of this section is to show that for 
sub- Poisson pp (i.e, pp that are dcx-smaller than Poisson pp), the critical 
intensity for /c-percolation of the Boolean model is non-degenerate. The re- 
sult for k = 1 (i.e., the usual percolation) holds under a weaker assumption 
of ordering of void probabilities and factorial moment measures. 

Given a pp of germs $ and the distribution of the typical grain G, ex- 
tending the definition of the Boolean model, we define the coverage field 
V$(x) := X^XiG* -'-[^ & Xi + Gi], where given {Gi} is a sequence of i.i.d. 
racs, distributed as G. The fc-covered set is defined as Gk{^,G) := {x : 
V$(x) > k}. Note that Ci(<I>, G) = C{^, G) is the Boolean model considered 
in Section 4. For /c > 1, define the critical radius for k-percolation as 

r^($,G) := inf{r : P (Cfc($,GW) percolates) > 0} , 

where, as before, percolation means existence of an unbounded connected 
subset. Clearly, rl{^,G) = rc($,G) < r^($,G). A gain, we shall abbreviate 
Cfc($,i?o(r)) by Cfc(<I>,r) and rl{^,{0}) by r,^($). 

Proposition 5.3. Let ^ be a stationary pp. For k > 1, X > 0, there 
exist constants c(A) and c{X,k) (not depending on the distribution of ^) 
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P \^^V<s>{zi) > nhj 

e-^^'^E^expjs^F^O 
e-^'^'^E(expL^y<i,,( 



i=l 



-snh 



exp|-A / (l-e*^-i^(^'^'))dxl , 



such that < c(A) < rl{^) provided $ <idcx ^\ and r^{^) < c{X,k) < oo 
provided <I> <ddcx ^x- Consequently, for $ <dcx $A combining both the above 
statements, we have that 

< c(A) < r^($) < r^{^) < c(A, k) < oo. 

Remark 5.4. 1. More simply, the theorem gives an upper and lower 
bound for the critical radius of a sub-Poisson pp dependent only on its 
mean measure (as this determines the A in ^x) and not on the finer 
structure. 

2. The following extensions follow by obvious coupling arguments. For a 
racsG C BoiR) a.s. with > 0, we have that rc(4>, G) > rc($,So(i?)) > 
(c(A) - R)+ and for a racs G D BoiR) a.s. with R > 0, r^{^,G) < 
r'^{^,Bo{R)) < {c{X,k)-R)+. 

Proof of Proposition 5.3. In order to prove the first statement, let 
idcx '^x &T^d r > 0. Consider the close packed lattice L*'^(2r). Define 
the response function lr{x,y) := l[x £ Qr{y)] and the corresponding shot- 
noise field V£{z) on L*'^(2r). Note that if C($,r) percolates then Zi.{V^, > 
1) percolates as well. We shall now show that there exists a r > such 
that Z2^(VJ, > 1) does not percolate. For any n and Zi E rZ'^, I < i < n, 
J2i=i lr{x^ Zj) = 1 iff X € Ur=i Qri^i) and else 0. Thus, from Lemma 5.2, we 
have that 

P(Vi(zO > 1,1 < i < n) < e-""exp|Ay" (e'^-i - l)dx| , 

= e-^"exp|A||UQ,(z,)ll(e^-l)|dx, 

(15) = (exp{-(. + (l-e^)A(2r)'^)}r. 

Choosing s large enough that e"** < (S'^ — 1)^^ and then by continuity of 
(s + (1 — e'*)A(2r)'^) in r, we can choose a c(A) > such that for all r < c(A), 
exp{ — (s + (l — e'*)A(2r)'^))} < (3*^ — 1)"-*^. Now, using the standard argument 
involving the expected number of open paths (cf Remark 5.1), we can show 
non-percolation of Ij2r{V^, > 1) for r < c(A). Hence for all r < c(A), C($, r) 
does not percolate and so c(A) < rd^). 

For the second statement, let ^ <ddcx ^x. Consider the close packed 
lattice L*'^(^). Define the response function lj.{x,y) := \[x G (5_£_(y)] 

and the corresponding additive shot-noise field V^{z) on L*'^(^). Note that 

Cfc(<I>,r) percolates \i7A_{V^, > \k/2\) percolates, where [a] = minjz G Z : 
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z > a}. We shall now show that there exists a r < cxd such that Z'^(V^, > 

\k/2\) percolates. For any n and Zi,l < i < n, from Lemma 5.2, we have 
that 

PiViizi) < \k/2] - 1,1 < i < n) 



,sn(rA;/2l-l) 



exp\x\\\jQ_T_{zi)\\{e-' -l)dx 



i=l 



(16) = (exp{-((i - e-nK^r - siim - mr. 

For any s, there exists c{X,k,s) < oo such that for all r > c{X,k,s), the 
last term in the above equation is strictly less than (S'^ — 1)~". Thus one 
can use the standard argument involving the expected number of closed 
contours around the origin (cf Remark 5.1) to show that Z'^(V^, > [/i;/2]) 

percolates. Further defining c(A, k) := infs>o c(A, k, s), we have that Ck{^, r) 
percolates for all r > c(A, k). Thus r'^i^) < c(A, k). □ 

For k = 1; i.e., for the usual percolation in Boolean model, we can avoid 
the usage of exponential estimates of Lemma 5.2 and work with void prob- 
abilities and factorial moment measures only, as in Section 4. The gain is 
two-fold: we extend the result to weakly sub-Poisson pp (cf. Section 3.3) and 
moreover, improve the bounds on the critical radius. 



Proposition 5.5. Let ^ be a stationary pp of intensity A and v-weakly 
sub-Poisson (i.e., it has void probabilities smaller than those of ^x). Then 

rem < C{X) := Vd (M^) < c(A, 1) < oo. 



Proof. As in the second part of the proof of Theorem 5.3, consider 
the close packed lattice L*^(;^). Define the response function lr{x,y) := 
l[x £ (5_2_(y)] and the corresponding extremal shot-noise field U^{z) := 
suPxe* ^r{z, X) on L*'^(^). Now, note that C{^, r) percolates if {z : UJ^{z) > 
1} percolates on 'L*'^{-^). We shall now show that this holds true for r > 
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c(A). Using the ordering of void probabilities we have 
PiUUz^) = 0,l<i<n) = P 



(17) 




< 

■t=i 

Clearly, for r > c(A), the exponential term above is less than (3*^ — 2)~^ and 
thus {z : U^{z) > 1} percolates by Peierls argument (cf Remark 5.1). It is 
easy to see that for any s > 0, exp{— A(;^)'^} < exp{ — (1 — e~'')A(;^)'^} and 
hence c(A) < c(A, 1). □ 

Proposition 5.6. Let ^ be a stationary, a-weakly sub-Poisson pp of 
intensity A (i.e., it has all factorial moment measures smaller than those of 
'^x). Then r,{^) > Ij^^^^j^^ > 0. 

Proof. We shall use the same method as in the first part of Theorem 5.3 
but just that we will bound the level crossing probabilities by using the fac- 
torial moment measures. As in Theorem 5.3, consider the close packed lattice 
L*'^(2r), the response function lr{x,y) := l[x S Qr{y)] and the correspond- 
ing shot-noise field V^{z) on L*'^(2r). We know that C(<I>,r) percolates only 
if Z2j.{V^, > 1) percolates. Let us disprove the latter for r < | (A(3d-i))i/d • 

P{VUzi)>l,l<i<n) = P{<^{Qr{zi))>l,l<i<n) 



< 



\i=l 
/ n 



< E[ll^x{Qr{z^)) 



\i=l / 

= (A(2r)^)" 

We can see that for r < ^ (A(3''-i))i/d ^2r(^$' — 1) does not percolate 

(cf Remark 5.1). This disproves percolation in C(<I>,r). □ 

Corollary 5.7. Combining the results of Propositions 5.6 and 5.5 we 

have < ^ (A(3d-i))i/d — '"c(*^) < Vd (^ ^°^^\ < oo for all weakly sub- 

Poisson pp. Examples of such pp are determinantal pp with the trace-class 
integral kernels (see Section 6.3). 
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5.3. Percolation in SINR graphs. Study of percolation in the Boolean 
model C(<I>,r) was proposed in Gilbert [1961] to address the feasibility of 
multi-hop communications in large "ad-hock" networks, where full connec- 
tivity is typically hard to maintain. The Signal-to-interference-and-noise 
ratio (SINR) model (see Baccelli and Blaszczyszyn [2001], Dousse et al. 
[2006] is more adequate than the Boolean model in the context of wireless 
communication networks as it allows one to take into account the interfer- 
ence intrinsically related to wireless communications. For more motivation 
to study SINR model, refer Blaszczyszyn and Yogeshwaran [2010] and the 
references therein. 

We begin with a formal introduction of the SINR graph model. In this 
subsection, we shall work only in M^. The parameters of the model are 
non-negative numbers P(signal power), A^(environmental noise), 7, T(SINR 
threshold) and an attenuation function £ : — )• M_|_ satisfying the following 
assumptions: i{x,y) = l{\x — y\) for some continuous function / : M+ — )• 
M+, strictly decreasing on its support, with ^(0) > TN/P, /(•) < 1, and 
xl{x)dx < 00. These are exactly the assumptions made in Dousse et al. 
[2006] and we refer to this paper for a discussion on their validity. 

Given a pp <I>, the interference generated due to the pp at a location x 
is defined as the following shot-noise field /$(x) := ^xe*\{x}KI^ ~ ^\)- 
Define the SINR value as follows : 

(18) 5/iVi?(x,y,cI>,7):^ PK\x - y\) 



N + -fPI^\{^}{y)- 

Let and be two pp. Let P, iV, T > and 7 > 0. The SINR graph 
is defined as G(<I>b, 7) := {^b, E{^B,^in)) where E{^b,^i,7) ■= 
{{X,Y) E : SINR{Y,X,^j,-f) > T,SINR{X,Y,^i,-f) > T}. The 
SNR graph(i.e, the graph without interference, 7 = 0) is defined as G{^b) ■= 
{^b,E{<!>b)) where E{^b) := {{X,Y) G : SINR{X,Y,^b,0) > T}. 

Observe that the SNR graph G(<I>) is same as the graph C(^>,r/) with 
2ri = Also, when = 0, we shall omit it from the parameters of 

the SINR graph. Recall that percolation in the above graphs is existence of 
an infinite connected component in the graph-theoretic sense. 

5.3.1. Poissonian back-bone nodes. Firstly, we consider the case when 
the backbone nodes (<5_b) form a Poisson pp and in the next section, we 
shall relax this assumption. When ^b = the Poisson pp of intensity 



^^The name shot-noise germ-grain process was also suggested by D. Stoyan in his private 
communication to BB. 
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A, we shall use G{X, and G(A) to denote the SINR and SNR graphs 

respectively. Denote by Ac(r) := X[rc{^x)/r)'^ the critical intensity for per- 
colation of the Boolean model C($>,i The following result guarantees the 
existence of a 7 > such that for any sub-Poisson pp <^ = $7, G(A, $,7) 
will percolate provided G(A) percolates i.e, the SINR graph percolates for 
small interference values when the corresponding SNR graph percolates. 

Proposition 5.8. Let A > Ac(r;) and <I> <idcx for some /i > 0. Then 
there exists a 7 > such that G{\, $, 7) percolates. 

Note that we have not assumed the independence of ^ and $a- In partic- 
ular, $ could be $A U where $0 is an independent sub-Poisson pp. The 
case $0 = was proved in Dousse et al. [2006]. Our proof follows their idea 
of coupling the continuum model with a discrete model. As in Dousse et al. 
[2006], it is clear that for = 0, the above result holds with Xc{ri) = 0. 

Sketch of the proof of Proposition 5.8. Our proof follows the ar- 
guments given in Dousse et al. [2006] and here, we will only give a sketch of 
the proof. The details can be found in in [Yogeshwaran, 2010, Section 6.3.4]. 

Assuming A > Xcipi), one observes first that the graph G(A) also per- 
colates with any slightly larger constant noise N' = N + 6' , for some 
5' > 0. Essential to the proof of the result is to show that the level-set 
{x : I^j{x) < M} of the interference field percolates (contains an infinite 
connected component) for sufficiently large M. Suppose that it is true. Then 
taking 7 = 6' /M one has percolation of the level-set {y : 7/<i>j(y) < 6'}. The 
main difficulty consists in showing that G{\) with noise N' = N + 5' per- 
colates within an infinite connected component of {y : /<i)j(y) < 6'}- This 
was done in Dousse et al. [2006], by mapping both models G(A) and the 
level-set of the interference field to a discrete lattice and showing that both 
discrete approximations not only percolate but actually satisfy a stronger 
condition, related to the Peierls argument. We follow exactly the same steps 
and the only fact that we have to prove, regarding the interference, is that 
there exists a constant e < 1 such that for arbitrary n > 1 and arbitrary 
choice of locations xi, . . . ,Xn one has P {I^j (xj) > M, z = 1, . . . , n) < e". In 
this regard, we use the first statement of Lemma 5.2 to prove, exactly as 
in Dousse et al. [2006], that for sufficiently small s it is not larger than 
for some constant K which depends on A but not on M. This completes the 
proof. □ 

5.3.2. Non-Poissonian back-bone nodes. We shall now consider the case 
when the backbone nodes are formed by a sub-Poisson pp. In this case, 
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we can give a weaker result, namely that with an increased signal power 
(i.e, possibly much greater than the critical power), the SINR graph will 
percolate for small interference parameter 7 > 0. 

Proposition 5.9. Let ^ be a stationary, u-weakly sub-Poisson pp and 
l^idcx ^^J. for some fj, > and also assume that l{x) > for all x E M_|_. 
Then there exist P, 7 > such that G($,$/,7) percolates. 

As in Theorem 5.8, we have not assumed the independence of and ^. 
For example, = <I> U <I>o where $ and $0 are independent sub-Poisson pp. 
Let us also justify the assumption of unbounded support for /(.). Suppose 
that r = sup{x : l{x) > 0} < 00. Then if C{^,r) is sub-critical, G(<I>,$/,7) 
will be sub-critical for any $7,^,7. 

Sketch of the proof of Proposition 5.9. In this scenario, increased 
power is equivalent to increased radius in the Boolean model corresponding 
to SNR model. From this observation, it follows from Proposition 5.5 that 
with possibly increased power the associated SNR model percolates. Then, 
we use the approach from the proof of Proposition 5.8 to obtain a 7 > 
such that the SINR network percolates as well. The details can be found in 
[Yogeshwaran, 2010, Section 6.3.4]. □ 

For further discussion on dcx ordering in the context of communication 
networks see Blaszczyszyn and Yogeshwaran [2010]. 

6. Examples of ordered point processes. 

In this section, we will give examples of pp for which some formal comparison 
of clustering is possible. In particular, in Section 6.1, we will give examples 
of families of perturbed lattice pp which are proved to be monotone with 
respect to the dcx order. In Section 6.1.1, we will show numerical evidences 
supporting the conjecture that within this class of pp the critical radius 
Tc is monotone with respect to the dcx order on the underlying pp. It is 
known that these pp can be considered as toy models for determinantal and 
permanental pp, which we show in Section 6.3 to be, respectively, weakly 
sub- and super-Poisson. Poisson-Poisson cluster pp are known to be dcx 
larger than Poisson pp. An example of such a pp - with Tc = - will be 
given in Section 6.2. This invalidates the conjecture on the monotonicity of 
Tc with respect to the dcx order of pp, in full generality. 

6.1. Perturbed point processes. Let $ be a pp on M'^ and •), •) 
be two probability kernels from M"' to non- negative integers Z+ and M'^, 
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respectively. Consider the following independently marked version of the pp 
= {{X,Nx,Yx)}x(^'i> where given 

• Nx, X G ^ are independent, non- negative integer- valued random vari- 
ables with distribution P {Nx G • | = M{X, •), 

• Yx = {Yix : i = 1,2, . . X £ ^ are independent vectors of i.i.d. ele- 
ments of M'^, with Yix^s having the conditional distribution P {Yix S • | ^I') 
X{X,-), 

• the random elements Nx,Yx are independent for all X G 
Consider the following subset of M'^ 

Nx 

(19) $^^'^*= U \J{X + Y,x}, 

xe$i=i 

where the inner sum is interpreted as when Nx = 0. The set ^p^*"* can 
(and will) be considered as a pp on M'^ provided it is locally finite. In what 
follows, in accordance with our general assumption for this article, we will 
assume that the mean measure of <^J"^'"* is locally finite (Radon measure) 

(20) / n{x)X{x, A-x) a{dx) < oo, for ah bBs B C R'^, 

where q(-) = E(<I>(-)) stands for the mean measure of the pp <I> and n{x) = 
YlT=i ^^{x, {k}) is the mean value of the distribution J\f[x, •). 

The pp (^P^'^^ can be seen as independently replicating and translating 
points from the pp with the number of replications of the point X G 
<I> having distribution M{X, •) and the independent translations of these 
replicas from X by vectors having distribution X{X, •). For this reason, we 
call a perturbation of <I> driven by the replication kernel N and the 

translation kernel X . 

An important observation for us is that the operation of perturbation of $ 
is dcx monotone with respect to the replication kernel in the following sense. 

Proposition 6.1. Consider a pp ^ with Radon mean measure a(-) 
and its two perturbations <I>^^'^* j = 1,2 satisfying condition (20), having 
the same translation kernel X and possibly different replication kernels Mj, 
j = 1,2, respectively. If M{x,-) <cx AA(x,-) (convex ordering of the condi- 
tional distributions of the number of replicas) for a-almost all x G M.^, then 

<dc. $r*- 

Proof. We will consider some particular coupling of the two perturba- 
tions <^f\ j = 1, 2. Given ^> and Yx = (Yix : i = 1, . . .) for each X G 
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A''-' 

let ^jx = [jiJi^{X + Yix}, where N^^ distribution J\fj{X, ■), j = 1,2, re- 



spectively. Thus <I>^'^'^* = X^xe* j — ^i"^ two considered pertur- 

bations. Note that given ^^"^ can be seen as independent superpositions 
of ^jx for X G Hence, by [Blaszczyszyn and Yogeshwaran, 2009, Propo- 
sition 3.2(4)] (superposition preserves dcx order) and [Miiller and Stoyan, 
2002, Theorem 3.12.8] (weak and Li convergence jointly preserve dcx order), 
it is enough to show that conditioned on ^ix ^dcx ^2X for every X S $. 
In this regard, given #, consider X £ ^ and let Bi,...,Bk be mutually 
disjoint bounded Borel subsets and / : R'^ — )• R, a dcx function. Define a 
real valued function (7 : Z — )• R, as 



g{n) := E /(sgn(n) J^(l[y.x G Si - X], ... , l[Y,x e Bt - X]) 



where sgn(?7-) = for n 7^ and sgn(O) = 0. By Lemma A.l, g{-) is a convex 
function on Z and by Lemma A. 2 it can be extended to a convex function g{-) 
on R. Moreover, E(5(iV^)|$) = E(^g{N'^)\<^) = E(/($,x(i?i), • • • , ^>jx(i?fc))|^ 
for j = 1, 2. Thus, the result follows from the assumption Nj^ <cx X'j^. □ 

Remark 6.2. The above proof remains valid for an extension of the 
perturbation model in which the distribution X{X, •) of the translations Yix 
depends not only on the location of the point X G $ but also on the entire 
configuration X{X,-) = Af(X, <!>,•), provided condition (20) is replaced 
by finiteness of Jj^jd J^in{x)X{x,4>,A — x) C{d{x,(j))), where C{d{x,(j))) is 
the Campbell measure of 

Perturbed pp, by Proposition 6.1, provide many examples of pp compa- 
rable in dcx order. We will be particularly interested in the following two 
special cases for the choice of the initial pp 

Example 6.3. Perturbation of a Poisson pp. Let $ be a (possibly in- 
homogeneous) Poisson pp of mean measure a{dx) on R'^. Let Af{x, •) = £! = 
1(1 G •) be the Dirac measure on Z"*" concentrated at 1 for all x G R'^ and 
assume an arbitrary translation kernel X satisfying a^^'^^{A) = J^a X{x, 
A — x) a{dx) < 00 for all bBs A. Then by the displacement theorem for 
Poisson pp, ^P^^'^ is also a Poisson pp with mean measure a^*^^*(dx). Assume 
any replication kernel N2{x, ■), with mean number of replications n2{x) = 
X^fcLi kN2{x, {k}) = 1 for all x G R"^. Then, by the Jensen's inequality and 
Proposition 6.1, one obtains a super-Poisson pp ^^^^ . In the special case, 
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when A/'2(x, •) is the Poisson distribution with mean 1 for all x gM.^, ^^^^ is 
a Poisson-Poisson cluster pp which is a special case of a Cox (doubly stochas- 
tic Poisson) pp with (random) intensity measure A(A) = Ylxi^^ X{x^A — x). 
The fact that it is super-Poisson was already observed in Blaszczyszyn and 
Yogeshwaran [2009]. Note that for a general distribution of 4?, its perturba- 
tion $2*^'^* is also a Cox pp of the intensity A given above. Other Poisson- 
Poisson cluster processes, with J\f2{x,-) not necessarily of mean 1, will be 
considered in Section 6.2. 

Example 6.4. Perturbation of a deterministic lattice. Assuming a de- 
terministic lattice <I> (e.g. $ = X'^) gives rise to the perturbed lattice pp 
of the type considered in Sodin and Tsirelson [2004]. Surprisingly enough, 
starting from such a one can also construct a Poisson pp and both super- 
and sub-Poisson perturbed pp. In this regard, assume for simplicity that 
$ = Z"^, and the translation kernel X{x, •) is uniform on the unit cube 
[0, 1)*^. Let Af{x, •) be the Poisson distribution with mean A {Poi{X)). It is 
easy to see that such a perturbation <I)P'^''* of the lattice Z*^ gives rise to a 
homogeneous Poisson pp with intensity A. 

Sub-Poisson perturbed lattices. Assuming for A/i some distribution convexly 
(ex) smaller than Poi{X), one obtains a sub-Poisson perturbed lattice 
pp. Examples are hyper- geometric HGeo{n, m, k), m,k < n, km/n = A 
and binomial Bin{n, X/n), X < n distributions which can be ordered 
as follows: 

HGeo{n,m, Xn/m) <cx Bin{m, X/m) <cx Bin{r, X/r) <cx Poi{X) 

for A < m < min(r, n); cf. Whitt [1985]^^. Specifically, taking Mi{x., •) 
to be Binomial Bin{n, X/n) for n > A, one obtains a dcx monotone 
increasing family of sub-Poisson pp. Taking A = n = 1 (equivalent to 
N'{x, •) = El), one obtains a simple perturbed lattice that is dcx smaller 
than the Poisson pp of intensity 1. A sample realization of this latter 
process (with $ being the unit hexagonal lattice on the plane rather 
than the square lattice) is shown in Figure 1. 
Super-Poisson perturbed lattices. Assuming for A/2 some distribution con- 
vexly larger than Poi{X), one obtains a super-Poisson perturbed lat- 

^'^Bin{n,p) has probability mass function PBin(n,p){i) = (")p'(l {i — 0, . . . ,n). 

HGeo{n,m,k) has probability mass function PHGeo(n.m,k){i) = (T) ("fc-T) / (fc) (max(fc - 
n + m,0) < i < m). 

^''One shows the logarithmic concavity of the ratio of the respective probability mass 
functions, which implies increasing convex order and, consequently, cx provided the dis- 
tributions have the same means. 
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tice. Examples are negative binomial NBin{r, p) distribution with rp/ (1 — 
p) = A, geometric Geo{p) distribution with 1/p — 1 = A; with 

Poi{\) <cx NBin{r2, A/(r2 + A)) <cx NBin{ri, X/{n + A)) 

<cx Geo(l/(l + A)) Geoipj) 

j 

with ri < r2, < Aj < 1, Xj = 1 and ^j/Pj = A + 1, where the 
largest distribution above is a mixture of geometric distributions hav- 
ing mean A; cf. Whitt [1985]. Specifically, taking A/'2(a;, •) to be negative 
binomial NBin{n, A/(n+A)) for n = 1, . . . one obtains a dcx monotone 
decreasing family of super- Poisson pp. Recall that NBin{r,p) is a mix- 
ture of Poi{x) with parameter x distributed as a gamma distribution 
with scale parameter p/{i — p) and shape parameter r. 

Remark 6.5. From [Meester and Shantliikumar, 1993, Lemma 2.18], we 
know that any mixture of Poisson distributions having mean A is cx larger 
than Poi{X). Thus, the super-Poisson perturbed lattice with such a repli- 
cation kernel (translation kernel being the uniform distribution on the unit 
cube) again gives rise to a Cox pp. A special case of such a Cox pp with M2 
being a mixture of two Poisson distributions was considered in Blaszczyszyn 
and Yogeshwaran [2009] (and called Ising-Poisson cluster pp). However, the 
proof of the fact that it is super-Poisson was based on the observation that 
the (random) density of this Cox pp is a conditionally increasing field. This 
argument can be extended to the case when the replication marks Nx, 
X G ^ constitute a field that is 1-monotonic ( [Grimmett, 2006, Ch. 2]). 
Due to space constraints, we refer to [Yogeshwaran, 2010, Ch. 5] for the 
details. A sample realization of the Cox pp obtained by the analogous per- 
turbation of the hexagonal lattice on the plane with A/2 being Bernoulli 
Bin{l, 1/5) mixture of Poi{5i) for i = 0, 1 is shown on Figure 1. 

Our interest in sub-Poisson perturbed lattices stems from their relations 
to zeros of Gaussian analytic functions (GAFs) (see Peres and Virag [2005], 
Sodin and Tsirelson [2004]). More precisely, Sodin and Tsirelson [2006] shows 
that zeros of GAFs have the same distribution as the pp UzeZ'^i-^ + -^z} 
for a Z^-shift invariant sequence {Xz}z^z'^- Simulations and second-moment 
properties (Peres and Virag [2005]) indicate that the zero set of GAFs exhibit 
less clustering (more spatial homogeneity) than a Poisson pp. The above 
example when seen in the light of the above-mentioned papers, asks the 

'^PG.o(r,){i)=p\l-pr-\PNB,n(r,Mi) + i^-PY- 
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0.50 0.55 0.60 0.65 0.70 

radius r 

Figure 2. Mean fractions of nodes in the two largest components of the sub-Poisson 
Boolean model C{$^[^{n),r) as a function of r; see Section 6.1.1. The underlying pp 
$^'^'^(n) are dcx increasing in n to the Poisson pp $a of intensity X = 2/(\/3) = 1.154701. 
The right dashed vertical line corresponds to the radius r = 0.5576495 which is believed 
to be close to the critical radius rc{^\) for the Poisson pp. The left dashed vertical line 
correspond to the critical radius rd^) — 0.5 for the non-perturbed lattice. 

question that we are currently not able to answer, whether zeros of GAF 
are sub-Poisson. 

We have seen above two examples of Cox pp (being some perturbations 
of a Poisson pp or a lattice) that are super-Poisson. A general class of Cox 
pp called the Levy-based Cox pp were proved to be super-Poisson in Bla- 
szczyszyn and Yogeshwaran [2009]. 

6.1.1. Numerical comparison of percolation for perturbed lattices. Con- 
sider Boolean model C(^>,r) with spherical grains of radius r. Recall that 
the critical radius Tc = rc{^) is the smallest radius r such that C(<I>, r) perco- 
lates with positive probability. In what follows, we will show some example 
of simulation results supporting the hypothesis that $1 <dcx ^2 implies 
rd^i) < rc(^>2). 

Consider two families of dcx ordered pp constructed as perturba- 

tions of the hexagonal lattice $ (of inter- node distance 1) with node trans- 
lation kernel X(x,-) being uniform within the hexagonal cell of x G $ and 
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n = 5 



n = 6 



n = 7 



n = 8 



Figure 3. Illustration of the phase transition in the clustering parameter n for the sub- 
Poisson Boolean model C {^^[J^ (n) , r) with r = 0.55; see Section 6.1.1. The largest com- 
ponent in the simulation window is highlighted. Bar-plots show the fraction of nodes m ten 
largest components. 




— I T"^ \ 1 1 — 

0.50 0.55 0.60 0.65 0.70 

radius r 

Figure 4. Mean fractions of nodes in two largest components of super-Poisson Boolean 
model C {^^gl„{n) , r) in function ofr; see Section 6.1.1. The underlying pp 'I>]vbL("-) '^''^ 
dcx decreasing in n to Poisson pp $a of intensity A = 2/(\/3) = 1.154701. 

different replication kernels N{x, •) (see Example 6.4). Specifically, assume 
binomial Bin{n,l/n) and negative binomial NBin{n, 1 /{I -\-n)) distribution 
for J\f{x, •) with n > 1. Recall from Example 6.4 that the former assumption 
leads to dcx increasing in n family of sub- Poisson pp ^p*^*"* = ^^^^^{n) con- 
verging to Poisson pp (of intensity A = 2/{\/3) = 1.154701) when n — )■ oo, 
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0.50 0.55 0.60 0.65 0.70 

radius r 



Figure 5. From sub- to super-Poisson. Figures 2 and ^ together. 



while the latter assumption leads to dcx decreasing family of super-Poisson 
pp $pe''* = ^^Bini''^) converging in n to the same Poisson pp. The crit- 
ical radius rc{^\) for this Poisson pp is known to be close to the value 
r = 0.5576495; 

In order to get an idea about the critical radius, we have simulated 300 
realizations of the Boolean model C(<1>^^''*, r) for r varying from r = 0.5 to 
r = 0.7 in the square window [0, 50] The fraction of nodes in the two largest 
components in the window was calculated for each realization of the model 
for each r and the obtained results were averaged over 300 realizations of the 
model. The resulting mean fractions of nodes in the two largest components 

^^Two dimensional Boolean model with fixed grains of radius r = 0.5576495 and Pois- 
son pp of germs of intensity A = 2/(^3) = 1.154701 has volume fraction 1 — e"^'''' — 
0.6763476, which is given in Quintanilla and Ziff [2007] as an estimator of the critical value 
for the percolation of the Boolean model. See also bound given in Balister et al. [2005]. 
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as a function of r are plotted on Figures 2 and 4 for binomial (sub-Poisson) 
and negative binomial (super-Poisson) pp, respectively. The two families 
are compared in Figure 5. The obtained curves support the hypothesis that 
the clustering of the pp of germs negatively impacts the percolation of the 
corresponding Boolean models. Moreover, Figure 3 illustrates what could be 
called the "phase transition for percolation in the clustering parameter" n 
for $Z*(n). 

6.2. Super-Poisson point process with a trivial percolation phase transi- 
tion. The objective of this section is to show examples of highly clustered 
and well percolating pp. More precisely we show examples of Poisson-Poisson 
cluster pp of arbitrarily small intensity, which are super-Poisson, and which 
percolate for arbitrarily small radii. 

Example 6.6. [Poisson-Poisson cluster pp with annular clusters] Let 
$Q, be the Poisson pp of intensity a on the plane M^; we call it the process 
of cluster centers. Consider its perturbation ^a'^''^ (cf. Example 6.3) with 
the following translation and replication kernels. Let X{x,-) be the uniform 
distribution on the annulus Bo{R) \ Bo{R — S) centered at x of inner and 
outer radii R — 5 and R respectively, for some < 6 < R < oo; see Fig- 
ure 6. Let Af{x, •) be the Poisson distribution Poi{fi). The process IS 
a Poisson-Poisson cluster pp; i.e., a Cox pp with the random intensity mea- 
sure A(-) := fi X^xe*c« '^(^' ■ ~^)- By [Blaszczyszyn and Yogeshwaran, 2009, 
Proposition 5.2], it is a super-Poisson pp. More precisely, <dcx ^a'^'^, 
where $a is homogeneous Poisson pp of intensity A = a^. 

For a given arbitrarily large intensity A < oo, taking sufficiently small 
a,R, 6 = R and sufficiently large /u, it is straightforward to construct a 
Poisson-Poisson cluster pp ^q'^''^ with spherical clusters, which has an ar- 
bitrarily large critical radius rc for percolation. It is less evident that one 
can construct a Poisson-Poisson cluster pp that always percolates, i.e., with 
degenerate critical radius Tc = 0. 

Proposition 6.7. Let ^a'^'^ be a Poisson-Poisson cluster pp with an- 
nular clusters on the plane as in Example 6.6. Given arbitrarily small 
a,r > 0, there exist constants a,fi,6,R such that < a,iJ,,5,R < oo, the 
intensity a/i of ^a'^'^ is equal to a and the critical radius for percolation 
rd^a'^'^) < r. Moreover, for any a > there exists pp $ of intensity a, 
which is dcx -larger than the Poisson pp of intensity a, and which percolates 
for any r > 0; i.e., rd^) = 0. 
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Figure 6. Poisson-Poisson cluster process of annular cluster; cf. Example 6.6. 

Proof. Let a, r > be given. Assume 5 = r/2. We will show that there 
exist sufficiently large /i, -R such that rc{^a' < t where a = aj^i. In this 
regard, denote K := 2TiR/r and assume that R is chosen such that K is an 
integer. For a a > and any point (cluster center) Xi ^ let us partition 
the annular support Ax^{R^ 5) := Bxi{R)\Bxi{R—S) of the translation ker- 
nel Xi+X{Xi, ■) (support of the Poisson pp constituting the cluster centered 
at Xi) into K cells as shown in Figure 6. We will call Xj "open" if in each of 
the K cells of Ax^ {R, S) , there exists at least one replication of the point Xi 
among the Poisson Poi{^) (with a = a/n) number of total replications of the 
point Xi. Note that given each point Xi G $q, is open with probability 
p{R, /i) := (1 — e~'^/^)^, independently of other points of <I>q. Consequently, 
open points of form a Poisson pp of intensity ap{R, fi); call it $open- Note 
that the maximal distance between any two points in two neighbouring cells 
of the same cluster is not larger than 2(6 + 2'kR/K) = 2r. Similarly, the 
maximal distance between any two points in two non-disjoint cells of two 
different clusters is not larger than 2[5 + 2'kR/ K) = 2r. Consequently, if the 
Boolean model C{^openi ^o{R^^)) with annular grains percolates then the 
Boolean model C{^a'^'^,r) with spherical grains of radius r percolates as 
well. The former Boolean model percolates if and only if C{^open, Bo{R)) 
percolates. Hence, in order to guarantee rc{^a'^'^) < r, it is enough to chose 
R, fx such that the volume fraction 1 - e-"P(^'/^)'^^' = 1 - e-«'P(^.A')'^^VM ig 
larger than the critical volume fraction for the percolation of the spherical 
Boolean model on the plane. In what follows, we will show that by choosing 
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appropriate R, fi one can make p{R, fj,)R'^ / fi arbitrarily large. Indeed, take 



^ := ij,{R) = log = log i? - - log logR) . 

Then, as i? — )• oo 

p{R,Ix)RVii = ^(i_e-W(2^R))2^i?/r 
A* 

Rr y/[SgR\^^R/^ 



27r(logi?- iloglogi?) V R 

^0(l)+log iJ-log(27r(log R- 1 log log i?))-0{l) v^I^ 



OO . 



This completes the proof of the first statement. 

In order to prove the second statement, for a given a > 0, denote a„ : = 
0/2" and let r„ = 1/n. Consider a sequence of independent (super-Poisson) 
Poisson-Poisson cluster pp <!>„ = ^^^'^"'i^" with intensities A„ := anfJ^n = On, 
satisfying rd^n) < ''n- The existence of such pp was shown in the first part 
of the proof. By the fact that are super-Poisson for all n > and by [Bla- 
szczyszyn and Yogeshwaran, 2009, Proposition 3.2(4)] the superposition <I> = 
UJ^i is dcx-larger than Poisson pp of intensity a. Obviously rc(<I>) = 0. 
This completes the proof of the second statement. □ 



Remark 6.8. By Proposition 6.7 we know that there exists pp ^ with 
intensity a > 0, such that r^^) = and ^dcx where $a is homoge- 
neous Poisson pp. Since one knows that rc{^a) > so (I* is a counterexample 
for the monotonicity of Tc in dcx ordering of pp. 

6.3. Determinantal and permanental point processes. In this section, we 
will give some results related to the dcx comparison of determinantal and 
permanental pp with respect to the Poisson pp. Refer to Ben Hough et al. 
[2006] for a quick introduction to these pp and for a more elaborate reading, 
see Ben Hough et al. [2009]. To this end, we will now recall a general frame- 
work (see [Ben Hough et al., 2009, Chapter 4]) which allows us to study 
ordering of determinantal and permanental pp more explicitly. 

6.3.1. Integral kernels. Let K : M.'^ x M"^ — )- C (where C are complex 
numbers) be a locally square-integrable kernel, with respect to on M."^^ 
Then K defines an associated integral operator /C/j on L^(D, ^) as /Cd/(x) = 

^®i.e., \K{x,y)\^ fi{dx) fi{dy) < 00 for every compact D C R'* 
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f^K{x,y)f{y)fi{dy) for complex- valued, square-integrable / on Z) (/ G 
L^{D,n)). This operator is compact and hence its spectrum is discrete. 
The only possible accumulation point is and every non-zero eigenvalue 
has finite multiplicity. Assume moreover that for each compact D the op- 
erator /C^) is Hermitian ^'^ , positive semi-definite and trace-class; i.e., 
I'^i^l ^ where Aj^ denote the eigenvalues of JCd- By the positive 
semi-definiteness of /Cd, these eigenvalues are non- negative. Further, one 
can show (cf. [Ben Hough et al., 2009, Lemma 4.2.2]) that for each compact 
D, there exists a "version" Kd{x, y) of the kernel K, defined on D' C D such 
that fi{D \ D') = 0, having the same associated operator ICd on L?'{D,^) 
which is Hermitian and positive semi-definite.^'^ Specifically one can take 
K£,{x,y) = (t>f {x)4>^ {y) where 4>f{-) are the corresponding normal- 

ized eigenfunctions oi K-d- 

Example 6.9. Determinantal pp. A simple pp on is said to be a 
determinantal pp with a kernel K{x, y) with respect to a Radon measure ^ 
on M'^ if the joint intensities of the pp with respect to the fcth-order product 
of ^ satisfy . . . , x^) = det {^K{xi,Xj)^^^^ .^^ for all k, where 

i^'^j) Ki j<k stands for a matrix with entries and det ( • ) denotes the 
determinant of the matrix. Note that the mean measure of the determinan- 
tal pp (if it exists) is equal to a(-) = J K{x, x) fj.{dx). Clearly, one needs 
assumptions on the kernel K(x, y) for the above equation to define the joint 
intensities of a pp. In what follows, we shall assume that the kernel K is an 
integral kernel satisfying the assumptions of Section 6.3.1. Then, there exists 
a unique pp on W^, such that for each compact D, the restriction of 
is a determinantal pp with kernel Kd if and only if the eigenvalues 
oi JCd are in [0, 1]. This latter condition is equivalent to Xj' € [0, 1] for all 
compact D; cf. [Ben Hough et al., 2009, Theorem 4.5.5]. We will call this pp 
determinantal pp with the trace-class integral kernel K{x,y). 

Example 6.10. Permanental pp. Similar to determinantal pp, one says 
that a simple pp is a permanental pp with a kernel K{x, y) with respect to 
a Radon measure /x on if the joint intensities of the pp with respect to 
11®^ satisfy p'^^\xi, . . . ,Xk) = per(iC(a;j,Xj))-^^^^.^^ for all k, where per( • ) 

^^i-e-, /p f{x) ICDg{x) fijdx) = g{x)K:Df{x) fi{dx) for all f,g € h'^{D,fi) 
i«i.e., J^f{x)ICDfix)^^{dx) > 

^^i.e. KD{x,y) = K{x,y) for almost all x,y £ D 

■^''Recall, a kernel K{x, y) is Hermitian if K{x, y) = K{y, x) for all x,y € R'', where z is 
the complex conjugate of 2: G C. It is positive semi-definite X^*Lj Zi ^i^ij ^ 

for all Zi £ C, i — 1, . . . , k, k > 1. 
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stands for the permanent of a matrix. Note that the mean measure of the 
permanental pp is also equal to a(-) = J K{x,x) fi{dx). Again, will assume 
that K(x,y) is an integral kernel as in Section 6.3.1. Then, there exists a 
unique pp i^P^"^^ on W^, such that for each compact D, the restriction of 
^perm D is a permanental pp with kernel K£); cf. [Ben Hough et al., 
2009, Corollary 4.9.9]. We will call this pp permanental pp with the trace- 
class integral kernel K{x,y). From [Ben Hough et al., 2006, Proposition. 35 
and Remark 36], we also know that ^p^"^"^ is a Cox pp with intensity field 
|-F|^, where F is some complex Gaussian process on W^. 

The following properties hold true for determinantal and permanental pp 
with a trace-class integral kernel K{x,y). 

• <1>'^*^* is a-weakly sub-Poisson, while (J)?"^''™- is a-weakly super-Poisson; 
both comparable with respect to the Poisson pp with mean measure 
a(-) given by a{D) = fj^ K£){x, x) fi{d) = X^jAj^, where the sum- 
mation is taken over all the eigenvalues Aj^ oi JCd- The proof fol- 
lows by the fact that Kjj^x, y) is Hermitian and positive semi-definite. 
Hence, by the Hadamard's inequality, < det [Kj:){xi,Xj))^^^ .^^ < 

Y[i=i ^D{xi, Xi) which implies (3). For ^p'^^^^ the proof follows from 
the permanent analogue of the Hadamard's inequality (see Marvin 
[1964]). In fact, for the above comparison to be true it suffices for the 
kernel K{x, y) to be Hermitian and positive semi-definite. 

• It is known that for each compact D, (^'^^^{D) has the same distribu- 
tion as the sum of independent Bernoulli random variables Bin{l, \^) 
cf. [Ben Hough et al, 2009, Theorem 4.5.3], while $p^™(D) has the 
same distribution as the sum of independent geometric random vari- 
ables Geo(l/(l + Aj^)), where the summation is taken over all eigenval- 
ues Xf of Kd; cf. [Ben Hough et al., 2009, Theorems 4.5.3 and 4.9.4]. 
Consequently 

(21) ¥^\D) PoiiY, Af ) <cx $P^™(I)) , 

i 

with the left inequality holding provided <I>'^'^* exists (i.e.; A^ G [0, 1] 
for all compact D). 

• Inequalities (21) imply that <I>°''^* is v-weakly sub-Poisson, while ^P"^*"™ 
is u-weakly super-Poisson. Consequently, is weakly sub-Poisson, 
while ^P*^^"^ is weakly super-Poisson. 

In Proposition 6.11, we will strengthen (21) proving dcx ordering of finite- 
dimensional distributions of ^'^^^ and ^P'^^"^ on mutually disjoint simulta- 
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neously observable sets Di, . . . ,-Dfc. Simultaneous observability means that 
the eigenfunctions of /Cy . , restricted to Dj are also eigenfunctions of /C 
for every i = 1, . . . ,k. 

Proposition 6.11. Let^'^'^* and^P^"^ be, respectively, the determinantal 
and permanental pp with a trace-class integral kernel K and with ^'^^^ being 
defined only if the spectrum of )C^d is in [0, 1] . Denote by <I>^°* the Poisson 
pp of mean measure a{-) given by a{D) = Xj' for all compact D, where 
the summation is taken over all eigenvalues \^ of K,d- Let Di, . . . , be 
mutually disjoint, simultaneously observable (with respect to the kernel K) 
compact subsets ofW^ and D = |J Dj. Then 

Proof. From [Ben Hough et al., 2009, Prop. 4.5.9], we know that 

J 

i=i 

with J denoting the number of eigenvalues of /C/j (J = oo and allowed 
and in the latter case the sum is understood as 0) and j = 1, . . . , J 
are are independent vectors, with components $,j = (^j^i, . . . , ^j^fc) whose 
distributions are P {^j^i = for all / ^ = 1) = ^f^i and P{^j,i = for 
alH = I,. . .,k) = 1 — A^, where {Xj,i}j=i,...,j are the eigenvalues of ICd^ 
and Xj' = X^JLi Aj^j, j = 1, . . . , J are the eigenvalues of JCd- Due to the 
independence of ^j's and the assumption Ylj=i^j < ^xd (local trace-class 
property of JCd), it is enough to prove for each j, that $j <dcx Cj where 
Cj = {Cj,i, ■ ■ ■ ,Cj,k) is the vector of independent Poisson random variables 
Cj^i ~ Poi{Xj^i). In this regard, note that the random vectors and Cj can 
be coupled as follows: $j = X^ili^i; Cj = ^fLi^i where N ~ _Bm(l,Aj^), 
M ~ Poi{Xf) and given N,M, ^i, i > 1 are independent multinomial 
vectors Mn/(1, Xf^^/Xf, Xf^k/^f) Then, <dcx Cj follows from the 
fact that Bin(l,p) <cx Poi{p) (see Example 6.4) and Lemmas A.l, A. 2. 
This completes the proof of the inequality for the determinantal pp. 

Mul{k,pi, . . . ,pk) with < Pi < 1, J2i=iPi — li probability mass function 
PMui{n,pi,...,Pk){ni, ■ ■ . ,nfe) = • • •Pfe'" forni + . . . + Tifc = n and otherwise. 
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Regarding the permanental pp, we know from [Ben Hough et al., 2009, 
Theorem 4.9.7] that the distribution of 

J 

($P-(Z)i),...,$P-(D,)) ^^r,^, 

i=i 

where 77^-, j = 1,...,J are are independent vectors, with rfj = X^^^^'j, 
where K ~ Geo(l/(l + A^) and, given K, independent multinomial vectors 
i = 1, . . . as above. Similarly, as for the determinantal pp, the required 
inequality follows from the ordering, C,j <dcx Vj for h which follows 
from the fact that Poi{\) <cx Geo(l/(l + A) (see Example 6.4) and Lem- 
mas A.l, A. 2. This completes the proof. □ 

Remark 6.12. Note that the key observation used in the above proof 
is that the number of points in disjoint, simultaneously observable sets 
Di, . . . , Dk can be represented as a sum of independent vectors, which them- 
selves are binomial (for determinantal) or Poisson (for Poisson) or geometric 
(for permanental) sums of some further independent vectors. This is exactly 
the same representation as for the perturbed pp of Section 6.1 (available for 
any disjoint sets Si, ... , -B/t); cf the proof of Proposition 6.1. In both cases, 
this representation and Lemmas A.l, A. 2 allow us to conclude dcx ordering 
of the corresponding vectors. 

Example 6.13. Ginibre process. Let be the determinantal pp on M? 
with kernel K{{xi,X2), (^1,^2)) = exp[{xiyi+X2y2)+iix2yi-xiy2)], Xj,yj G 
M, j = 1,2, with respect to the measure iJ,{d{xi, X2)) = 7r~^exp[— — 
X2] dxidx2- This process is known as the infinite Ginibre pp. Denote by 
^''^ = {l-'^ip : Xi G ^^}, the pp on of the squared radii of the points of 

. We know that an arbitrary finite collection of the annuli centered at the 
origin Di = {{xi,X2) ■ Vi < x\ + x\ < Ri} is simultaneously observable for 
this pp; cf. [Ben Hough et al., 2009, Example 4.5.8]. Using this observation. 
Proposition 6.11 and the fact that dcx order of pp on M is generated by the 
semi-ring of intervals, we conclude that is dcx smaller than the Poisson 
pp of unit intensity on M+. A partial result, ^*^([0,r]) <cx $^([0,r]) for 
all r > 0, was proved in Blaszczyszyn and Yogeshwaran [2009]. 

7. Concluding remarks. We come back to the initial heuristic dis- 
cussed in the Introduction — clustering in a point process should increase 
the critical radius for the percolation of the corresponding continuum per- 
colation model. As we have seen, even a relatively strong tool such as the 
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dcx order falls short, when it comes to making a formal statement of this 
heuristic. 

The two natural questions are what would be a more suitable measure of 
clustering that can be used to affirm the heuristic and whether dcx order 
can satisfy a weaker version of the conjecture. 

As regards the first question, one might start by looking at other depen- 
dence orders such as super-modular, component-wise convex or convex order 
but it has been already shown that the super-modular and component- wise 
convex orders are not suited to comparison of clustering in point processes 
(cf. [Yogeshwaran, 2010, Section 4.4]). Properties of convex order on point 
processes are yet to be investigated fully and this research direction is in- 
teresting in its own right, apart from its relation to the above conjecture. In 
a similar vein, it is of potential interest to study other stochastic orders on 
point processes. 

On the second question, it is pertinent to note that sub-Poisson point 
processes surprisingly exhibited non-trivial phase transitions for percolation. 
Such well-behavedness of the sub-Poisson point processes makes us wonder 
if it is possible to prove a rephrased conjecture saying that any homogeneous 
sub-Poisson pp has a smaller critical radius for percolation than the Poisson 
pp of the same intensity. Such a conjecture matches well with [Benjamini 
and Stauffer, 2011, Conjecture 4.6]. 

APPENDIX A 

The following result, similar to [Meester and Shanthikumar, 1993, Lem- 
ma 2.17] is used in the proof of Propositions 6.1 and 6.11. 

Lemma A.l. Let = (^.^, . . . G M'^, = 1, 2, . . .) he independent, 
identically distributed vectors of (possibly dependent) non-negative random 
variables. Suppose f is a dcx function on R'^. Then, the function g defined 
on Z by g{n) = ^(^f{sgn{n)Y^\^^^^)^ for n ^ and g{0) = is convex 
on 7L. 

Proof. We will prove that g{n) has non- negative second differences 

(22) g(n - 1) + g(n + 1) - 2g{n) > for all n G Z 

and use the ffist part of Lemma A. 2. To prove (22), define G{n,m) := 
HZn+i for < n < m and G{n, n) := (0, 0, . . . , 0) G M'^ for n > 0. We 
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have for n > 1, 

2g{n) = 2E(/(G(0,n))) 



= E^/^G(0,n-l) + G(n-l,n)jj +E^/^G(0,n) 

= E(/(G(0,n-l) + G(n,n + l)) +E(/(G(0,n 

= E(/(G(0,n-l) + G(n,n + l) +/(G(0,n) 

< E(/(G(0,n-l)) +/(G(0,n) + G(n,n + l) 
= 9{n-l) + g{n + l) , 

where for the third equahty we have used mutual independence of G(0, n — 
1), G{n — 1, n), G(n, n + 1) and the fact that G{n — 1, n) and G(n, Ti + 1) have 
the same distribution, while the inequality follows from the dcx property of 
/ and the assumption > 0. This proves (22) for n > 1. Similar reasoning 
allows to show (22) for n < — 1. Finally, note that for n = 



2^(0) = 2/(^(0,...,0) 

= e(/(-G(0,1) + G(0,1)) +/((0,...,0) 

< e(/(-g(o,i))+/(g(o,i) 

= 5(-1) + 5(1), 



□ 



We will prove the following two technical results regarding convex func- 
tions. We were not able to find their proofs in the literature. 

Lemma A. 2. Let g{n) he a real valued function defined for all integer 
n ^"L and satisfying condition (22). Then the following two statements are 
true: 

1. For alln>2 

n n 

(23) ^ Y.^'9{ki) 

i=l 1=1 

for all ki £ Z and < Ai < 1, ^17=1 = 1 such that Yl'i=i ^i^i ^ ^■ 
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2. Function g{-) can be extended to a real valued convex function defined 
on real numbers M. 

Proof. Ad 1. As mentioned in [Marshall et al., 2009, Section V.lG.B.lO.a] 
it is easy to see that (22) is equivalent to (23) with n = 2. Assume now 
that (23) holds true for some n>2 (and all < Aj < 1, A;^ G Z, i = 1, . . . , n 
satisfying ^11=1 Aj = 1, Yll=i ^i^i ^ We will prove that it holds true for 
n + 1 as well. In this regard, define for a given /c € Z and distinct (otherwise 
we use directly the inductive assumption) /ci, . . . , G Z, the following 
functions: 



An+i — A„-|.i(Ai, . . . , A„_i) 
F(Ai, . . . , A„_i) 



kn ~ kn-\-l 

n-1 

:= 1 — Xj — An(Ai, . . . , Xn~l) 

1=1 

n-1 

:= ^ Xig{ki) + An(Ai, . . . , An-i)5'(^n) 
1=1 

+An+l(Al, . . . , An-l)5'(Wl) • 



Note that for any Ai, . . . , An-i we have Y^^=i Aj = 1 and XlILi ~ ^• 
Consider the following subset of the n — 1-dimensional unit cube 

C:= {(Ai,...,An_i) G [0,l]"-i :0< An< 1,0< An+i < l} . 

The proof of the inductive step will be completed if we show that F{-) > g{k) 
on C. In this regard note that C is closed and convex. Assume moreover that 
C is not empty; otherwise the condition (23) is trivially satisfied. Note also 
that F{-) is an affine, real valued function defined on M"~^. Hence, by the 
maximum principle, the affine (hence convex) function —F attains its max- 
imum relative to C on some point (A?, . . . , X^_i) G dC of the boundary of 
C. Consequently, we have F(-) > F{X\, . . . , A^_^) on C and the proof of the 
inductive step will be completed if we show that F{X\, . . . , A^_^) > g{k). In 
this regard, denote A° = An(A'j', . . . , A°_J and A°^_i = An+i(A?, . . . , A°_i). 
Using the continuity of the functions An(-) and An+i(-) is not difficult to ver- 
ify that (A?, . . • , X^_i) G dC implies A^ = for some j = 1, . . . , n + 1. Thus, 

by our inductive assumption, F{X^, . . . ,X^_i) = Yl^=i i^j "^^diki) > g{k), 
which completes the proof of (23) for all n > 2. 

Ad 2. We will recall the arguments used in Yan [1997] to show that a 
function satisfying (23) for all n > 2 (called globally convex function there) 
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has a convex extension on M. In this regard, consider the epigraph epi(5f) := 
{(fc, /i) G Z X M : /i > g{k)} of g and its convex envelope ep'f°{g). It is easy 
to see that epi'^°{g) = {{x,fi) G : /i > ^21=1 ^idi^i) some ki £ Z, 
< Aj < 1, X]r=i -^j = 1 and Ya=i ^i^i = Define ^(x) := inf{^ : {x, /i) G 
epi'^°(g')} for all x G M. The convexity of epi'^°{g) implies that g is convex 
on M and the global convexity (23) of g implies that g is an extension of g. 
This completes the proof. □ 
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